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ABSTRACT 


This manual describes the use of the computer program BUCLASP2. The 
program is intended for linear instability analyses of structures 
such as uni directionally stiffened panels. Any structure that has 
a constant cross section in one direction, that may be idealized as 
an assemblage of beam elements and laminated flat and curved plate 
strip elements can be analyzed. The loadings considered are 
combinations of axial compressive loads and in-plane transverse 
loads. The two parallel ends of the panel must be simply supported 
and arbitrary elastic boundary conditions may be imposed along any 
one or both external longitudinal side. 

This manual consists of instructions for use of the program with 
sample problems, including input and output information. The 
theoretical basis of BUCLASP2 and correlations of calculated results 
with known solutions, are presented in Reference 1, 
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1.0 SCOPE OF BUCLASP2 

The computer program BUCLASP2 is designed to perform instability analyses for 
structures such as uni di recti onal ly-sti ff ened , rectangular composite panels. 
The structure is idealized as an assemblage of beam elements laminated 
curved plate strip elements and laminated flat plate strip elements. Each 
element extends the full length of the structure. The element edges which 
are normal to the longitudinal axis of the panel are assumed to be simply 
supported and any external edge parallel to the longitudinal axis may be 
arbitrarily constrained by specifying the appropriate spring constants. 

Plate offsets can be employed when needed to correctly idealize a structure. 

The analysis is based on linear, elastic theory with prebuckling stresses 
and deformations ignored. Stress analyses cannot be performed by using 
BUCLASP2. The loadings considered are special combinations of constant in- 
plane axial and transverse loads and strains. 

The individual load dependent elemental stiffnesses are solved exactly 
and then appropriately merged to form the total stiffness matrix of the 
structure. Inter-element equilibrium and displacement compatibility 
conditions are automatically satisfied. By setting the determinant of 
the load dependent total stiffness matrix to zero, a nonlinear eigenvalue 
problem is created. An iterative scheme, based on successive guesses of 
the loads, is then used to extract the minimum eigenvalue. In order to 
keep the number of guesses to a minimum, the method described in Reference 
1 is used. 

The critical loads are solved for each specified axial wave number, m, and 
mode shape calculations are made only for the m value that corresponds 
to the lowest critical load. In order to get the lowest buckling load, 
the structure should be analyzed for different m values and a search 
made to find the lowest critical load of the structure. 


1.1 



2.0 BASIC INFORMATION 

This section presents the basic information required by BUCLASP2 users 
for performing buckling analyses. Figure 2.1 shows the cross section 
of an arbitrary structure. The circled numbers are element numbers. 
Extremities of each plate element and the shear centers of each beam 
element are indicated by dots. The dashed line is drawn through the 
mid-planes of each plate element. 
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2.1 Geometry 


2.1.1 Coordinate Systems 

All structural geometries are defined relative a global (overall) 
right-handed rectangular Cartesian coordinate system. The axes 
of this triad are denoted by Xg, Yg, and Zg in Figure 2.1. The 
longitudinal axis of the structure is in the same direction as X Q . 
Local element coordinate systems will be discussed individually 
for each element in Section 2.2. 

2.1.2 Nodes 

The dots shown in the Figure 2.1 will be referred to as nodes. They 
are actually lines that are parallel to Xg . A node is identified by 
its coordinates in the global system and by its user-assigned node 
number. This node number is only for the user's convenience since 
the program automatically numbers the nodes sequentially according 
to the input order. The first input node will be internal node 
number 1, the second input node, 2, etc. Generally, the user node 
numbers will not be the same as the internal node number. No 
reordering of nodes is done internally. 
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2.2 Elements 

There are three basic types of elements: (1) flat plate, (2) curved 

plate, and (3) beam elements. All elements extend the full length 
of the structure. A plate element is defined by a pair of node 
numbers that corresponds to the sides of the element. In Figure 2.2 
the plate element is defined by nodes A and B. A beam element is 
defined by a node that corresponds to the beam shear center. In Figure 
2.1 the following observations can be made: 

(a) © , ® , © . ® , ® and © are flat plate elements. 

(b) © , © and ® are curved plate elements. 

(c) © and © are beam elements. 
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2,2.1 Plate Elements 

There are two types of plate elements, flat and curved. The following 
discussion applies to both types except when reference is made to a 
particular type of plate element: 

Local Coordinate System 

The local coordinate system of a plate element is defined by the 
ordering of the node numbers that identify the element. The local 
y axis of the plate element is in the direction from the first node 
of the element to its second node. For example in Figure 2.2 the 
element is defined from node A to node B, and the resulting direction 
of the local y axis of the element is as shown. The local y axis of 
the curved plate element is along the arc from A to B. The local x 
axis is parallel to and in the same direction as X Q . The local z 
axis is determined by the right-hand rule. The origin of the local 
axis is on the mid-surface of the element and midway between its 
nodes. 

Properties 

In its most general form, the plate element can be multi-layered, each 
layer being orthotropic with respect to the local plate axes. In its 
simplest form, the element degenerates to a single isotropic layer. 

The thickness of each layer is assumed to be constant. The user must 
specify the thickness and material property of each layer. 

The first lamina of a laminated plate element is defined as the 
external layer in the negative z direction. The last lamina is 
the extreme layer in the positive z direction (see Figure 2.2). It 
is important that the layers are input in this particular order. 
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Radius of Curvature of Curved Plate 

The curved plate element is actually a strip of a circular 
cylindrical shell. Therefore, only one radius of curvature, 

R, need be defined for this element. The sign convention for 
R is defined in the following manner: R is positive if the 

vector from the element to the center of curvature is in the 
same direction as the local z axis of the element. The magnitude 
of R is measured from the mid-surface of the element to the center 
of curvature. In Figure 2.3 the sign of R for elements (a) and 
(b) is negative, whereas it is positive for element @. 
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2.2.2 Beam Elements 

There are two basic types of beam elements that are available in 
BUCLASP2. They are the laminated beam elements, circular and 
rectangular in cross section, and the general beam element. The 
local x axes of all beam elements are directed in the same direction 
as the global axis and parallel to it. 

Laminated Beams 
Rectangular: 

The rectangular laminated beam is shown in Figure 2.4(a). The 
depth of each layer is b. The origin of the local axis system 
is taken at the shear center, 0, of the beam. The local y axis 
of this element could either be in the direction as shown in Figure 
2.4(a) or in the opposite direction. The local z axis is determined 
by using the right-hand rule. 

The laminae are ordered sequentially in the direction of the local 
y axis as indicated by Figure 2.4(a). Thus, if the local y axis 
were in the opposite direction, the ordering of the layers would 
be reversed. 

Circular: 

The beam is laminated in the manner shown in Figure 2.4(b). The 
first layer is the innermost layer and the ordering of the layers 
extends outward. The local y axis of this element can be in any 
direction in the plane of the beam cross section. The local z axis 
is again determined by the right-hand rule. 
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(a) Laminated Rectangular Beam 


(b) Laminated Circular Beam 


Figure 2.4 Laminated Beam Elements 
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General Beam Element: 

The local y axis of the element can be in any direction in the 
plane of the beam cross section. The local z axis is defined by 
the right-hand rule. Properties such as inertias, areas, torsion 
constants and warping constants are input with respect to the 
local coordinate system. The orientation of a beam element is 
defined relative to the global coordinate system by specifying the 
angle between the global and local axes. 
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2.3 Plate Offsets 

The element stiffnesses are originally defined at the element "local 
nodes". A plate element has two "local nodes" defined as the longi- 
tudinal edges lying in the mid-surface and the "local node" for a beam 
element is defined by its shear center. When two or more elements are 
joined together, the appropriate element stiffnesses can be merged at 
the juncture if the "local nodes" of all elements are coincident at the 
juncture. If the "local nodes" at this juncture are not coincident, it 
is necessary to define a common node where all the element stiffnesses 
can be merged properly. By necessity this common node must be a user 
input node. By using plate offsets, the user effectively defines the 
element stiffness at some node other than the plate "local node", thereby, 
giving him a tool to meet the above requirements. This non-" local node" 
is referred to as an "offset node". It must, however, be one of the 
user-defined nodes. 

A plate with offsets can be interpreted as a pseudo plate element with a 
"link" connecting ene of the plate "local nodes" to an "offset node". 

This link allows the element stiffness at a "local node" to be transferred 
to an "offset node". 

There is no offset capability for the beam element in this program. This, 
however, is not detrimental since in all cases plate offsets can be used 
to idealize non-coincident junctures. 

A situation where plate offsets can be used is when the mid-surfaces of 
plates do not intersect at the juncture. Offsets can also be used when 
the shear center of a beam does not coincide with the mid-surface of an 
adjoining plate. 
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The plate offset is a vector which is specified by two components, y Q 
and z Q from the edge of the plate element (see Figure 2.5.1 and 2.5.2). 
The y Q and z Q vector components are measured relative to the local 
y and z axis respectively. They are positive when their directions 
coincide with the positive directions of the local axes. z Q is 
measured from the negative-z plate surface (exterior surface in the 
negative-z direction) and y Q is measured from the plate "local node". 

It must be emphasized that z Q is not measured from the "local node". 

Figure 2.5.1 illustrates the cross section of a panel in which two 
plate elements intersect a beam element. The input nodes coincide with 
the plate "local nodes" but they are not coincident with the beam shear 
center. Since these nodes are not coincident at the juncture and a 
beam offset capability is not available, offsets for plate elements (T) 
and (?) must be specified. Consequently, the element stiffnesses will 
be merged at 3, the beam node. 


Figure 2.5.2 illustrates two examples where plate elements of unequal 
thicknesses meet. In the top figure, the mid-surfaces of plate 
elements (T) and (3) with equal thicknesses passes through the input 
nodes. Element (2) , however, is a two-layered plate with a different 
thickness. Since the stiffness for element (T) is formulated relative 
to its "local nodes", offsets must be defined to transfer this stiffness 
to nodes 2 and 3 to ensure the desired structural integrity. If offsets 
were not specified for this example, the true configuration of this 
structure would not be retained in the analysis. That is, element (2) 
would be situated such that its mid-surface coincides with the mid- 
surfaces of (T) and (3). 

It should be noted that z q is not the offset distance from a plate 
"local node" to the "offset node". This offset component, however, is 
used by the program to formulate the actual offset- transfer of the 
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stiffness. For example, the stiffness at the mid-surface of ( 2 ) in 
Figure 2.5.2, is transferred to nodes 2 and 3 based on the z Q 
offset indicated therein. 
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©0 


y Q Positive ( for Element (T) ) 


y D Negative ( for Element (3) > 



and ( 3 ) Plate Elements 
Beam Element 


Figure 2.5.1 Offsets in Presence of Beams 





(T), ( 2 ) and ( 3 ) Plate Elements 

Figure 2.5.2 Plate Offsets 
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2.4 Boundary Conditions 

An external side (parallel to the X Q axis) of a plate element may 
be constrained by specifying the following boundary conditions; 

(a) Simply supported 

w = M^2 = ^22 = u = 0 

(b) Clamped 

w = e = n 22 = u = 0 

(c) Free 

1 ■ M 22 ' N 22 * N 12 * 0 

(d) Sprung 

The constrained or sprung degrees of freedoms are defined with 
respect to the local coordinate system of the plate element. 
Boundary conditions cannot be applied to a beam element node 
or an internal node. For a sprung node of a plate element, 
the spring stiffnesses for all four degrees of freedom must be 
specified. To delete or constrain a freedom, use a large value 
for the corresponding spring constant, e.g., 1 X 10^ 00 . 
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2.5 


Element Loadings 

The distribution of the total buckling load to individual elements 
is discussed in this section. In all cases the inplane loads 
and N 2 2 for plate elements and axial load for beam elements are 
uniform. Shear loads N-j 2 and moments are not considered. ''Affected" 
elements are herein defined as plate elements subjected to transverse 
N 22 load. The N 22 load is equal to zero for non-"affected" elements. 
All "affected" plate elements are subjected to the same N 22 load. 

Following are the special options that are available in the code. 

In the most general case, all options listed below can be used to 
solve buckling problems of biaxially loaded panels. For uniaxial 
buckling analysis, these options can be used by setting to zero, 
properly selected input parameters. 

(1) Biaxial loading, constant specified N 22 load 

This option is to be used for a panel that is loaded by a 
constant side load. The side load N 22 is applied to all 
"affected" elements. While preserving strain compatibility 
in the longitudinal direction, the critical load is 
solved. If the input N 22 load is specified as zero, the 
unixial buckling load of the panel will be computed. For 
a stiffened panel under biaxial loads, the stiffener 
webs and outstanding legs would normally not be subjected to the 
transverse N 22 load, while the skin of the panel would be loaded 
in the transverse direction. For this panel, the skin elements 
would be "affected" and the stiffener webs and outstanding legs 
would be non-"affected". 

(2) Biaxial loading, constant specified biaxial load ratio 

The biaxial load ratio is defined as the N n load of the first 
"affected" plate element divided by the N 22 side load. While 
preserving longitudinal strain compatibility of all elements, 
including beams, and by keeping the biaxial load ratio constant, 
the critical buckling load is calculated. 
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(3) Biaxial loading, constant specified longitudinal strain 

The specified longitudinal strain is the same for all elements in 
the panel and will remain constant, N 22 critical load will be 
computed for this option. Only "affected" plate elements will be 
subjected to the N 22 load. 

(4) Biaxial loading, constant specified axial load 

The panel buckling load will be established for a constant total 
axial load. The N-j-j load for the first plate element is specified. 
Initially, by assuming longitudinal strain compatibility at zero 
N 22 load, the axial loads in all elements are found. These 
axial loads will remain constant. The N 22 critical loads on 
"affected" elements will then be computed. At the buckling load, 
longitudinal strain compatibility will not necessarily be preserved. 
If N^-j is specified as zero, the panel buckling load due to side 
load only will be calculated. 

For the options listed above, all elements in the panel, beams 
and plates will be subjected to loads, the exception being the 
case where the specified N n load in option (4) is set equal to 
zero. Tensile transverse loads and strains can also be specified. 
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2.6 


Substructures 


Three types of substructures are considered in BUCLASP2/ (1) start, 

(2) repeat and (3) end. The bas i c reason for including these types 
of substructures is to efficiently analyze panels with many repeated 
structural groups along the panel width. Examples of these types of 
panels are shown in Figure 2.6. The merged stiffness matrix for this 
type of panel consists of many repetitive submatrices because the 
repetitive substructures yield identical stiffnesses. BUCLASP2 takes 
advantage of this repetitiveness by dividing the stiffness matrix into 
three distinct blocks , one for each start, repeat and end substructure, 
and uses only one repeat block (corresponding to the repeat substructure) 
for the buckling load computations. Thus the generation and solution 
times for the stiffness matrix are reduced and the storage requirements 
are also alleviated. 

The user may divide a panel, with or without repetitive structural 
groups, such that it consists of one start substructure, one, many, 
or no repeat substructures and one end substructure. Further, a 
panel may be idealized with the following combinations of substructures; 

(1) start substructure only 

(2) start and end substructures only 

(3) start, repeat (one or more) and end substructures. 

It should be obvious at this point that panels with no repeatable 
structural groups can have a maximum of one repeat substructure. 

It is recommended that the start substructure be larger of the 
three substructures, for efficiency reasons in the determinant 
calculations. 
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2.6.1 Flat Panel with Repeated Substructures 

For a panel with repeated substructures* a reduced version of the 
actual panel may be specified for input. For example, see Figure 2.7. 
(a) is the actual panel and (b) is the reduced panel used for input. 
The nodal coordinates of 7 to 10 of (b) are input such that they 

model the end substructure of the actual panel, (a), as if it is 

moved to the left and attached to the first repeat substructure. 
Therefore, nodes 7 to 10 of the reduced panel, (b), are related to 
nodes 23 to 26 of (a), the actual panel. Consequently, nodes 11 to 

22 of (a) are not directly included by coordinates in the input, but 

are implicitly specified through the repeat substructure. 

After drawing a sketch of the actual structure, the user must identify 
the nodes that belong in the three different substructure types. 

For the start and end substructures, the only limitation is that 
there must be at least one interior node, a node without any specified 
boundary conditions, in each substructure. To ensure that the repeat 
substructure is accurately identified, the substructure corresponding 
nodes must be identified. The corresponding nodes are those node 
pairs at which the stiffnesses are identical. In other words, all 
elements meeting at corresponding nodes in all repeat substructures 
must be identical in all respects, including loadings and boundary 
conditions. 

For the example from Figure 2.7 (a), inspection shows that nodes 7, 11, 
15 and 19 are the corresponding nodes to node 3. Therefore, node 3 
is a valid node for the repeat substructure. Similarly, node 6 is 
a valid node for the repeat substructure with nodes 10, 14, 18, 22 
as correspondence nodes; node 4 with 8, 12, 16, 20; and node 5 with 
9, 13, 17 and 21. One should note that node 23 is not a corresponding 
node of node 3 because the specified boundary condition at node 26 
results in the stiffness contribution of element 23-26 to node 23 
being unequal to the stiffness contribution of element 3-6 to node 3. 
Thus, node 23 may not be included in the repeat substructure. 
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From this user exercise one can conclude that, a) the substructure 
unit made up of nodes 3, 4, 5, 6 is a repeat substructure, and b) 
there are (5) repeat substructures. 

Node 2 is included in the start substructure because boundary 
conditions are specified for node 1. The selection of nodes for 
the end substructure is obvious at this point. 
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2.6.2 Curved Panel with Repeated Substructures 

For curved panels with repeatable substructures, repeatability of 
the nodes for the repeat substructure must be viewed in the curved 
panel cylindrical coordinate system which is shown In an example of 
Figure 2.8. R p , 0 and Xq are the axes of this coordinate system. 
The same guidelines as the ones discussed for flat panels are 
applicable for curved panels. The reduced structure for input 
is shown below the actual structure. As in flat panels, the end 
substructure is moved from its original position and is attached 
to the repeat substructure in the reduced version. The nodal 
coordinates are input in the global coordinate system (Xq, Yq and 
Z Q ) for the reduced structure, 

A0, the include angle, in Figure 2.8 is the angle measured from a 
point in a repeat substructure to its corresponding point in the 
next repeat substructure. The positive sign of A0 is as shown. 

If the curved panel in Figure 2,8 is to be repeated from right to 
left instead of left to right as shown, the sign of A0 would be 
negative. In terms of vectors, the positive rotation vector A0 
is parallel to and in the same direction as Xq using the right 
handed convention. 
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2.6.3 Repeat Substructure Interrelationship Node Pairs 

The mode shapes for elements of the panel (not the reduced 
structure) is provided by the program. These elemental mode 
shapes are derived from the relative displacements which were 
determined individually for each node of the actual panel 
(not just the reduced structure). The previous input information 
on repeat substructures was adequate to produce (internally) 
the total structure stiffness from the (input) reduced structure. 
However, additional information is required to identify elements 
that connect one repeat substructure to another repeat substructure 
for the purpose of calculating elemental mode shapes. 

From the definition of a repeat substructure, any element, that 
connects the start substructure to the adjoining repeat substructure, 
has a corresponding element (topologically similar), that connects 
the repeat substructure to the end substructure. The requirement 
then becomes one of identifying; 

(a) for each connecting element between the start and repeat 
substructures, the node that lies in the start substructure, 
and 

(b) for each connecting element (topologically corresponding to 
(a) ) between the repeat and end substructures, the node 
that lies in the repeat substructure. 

The nodes indicated above become logically a node pair between 
corresponding elements and are defined as "interrelationship 
node pairs". 

As an aid to the user, the following steps are recommended 
(see Figure 2.9). 

(1) After encircling the nodes for different substructures 
as shown in Figure 2.9, identify all element(s) that 
connects the start substructure to the repeat substructure. 
(Elements (j) , @ , Q) , © and @ of Figure 2.9) 
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(2) Identify all element(s) connecting the repeat and end 
substructures. (Elements ^3) , (0) , (0) , (0) and 
© of Figure 2.9) 

(3) Pair all elements of (1) above, with corresponding elements 

(topologically similar) from (2) above. For the panel in 
Figure 2.9, the element pairs are: ( ( 3 ) , © ) , 

( © , <©),(©. ©).(©. © ) and 

( © . © ). 

(4) For each corresponding element pair, select the start 
start substructure node from element(s) of (1) above 
and the repeat substructure node from element(s) of (2) 
above. These two nodes constitute the interrelationship 
node pair. For element pair ( (?) , © ), they are 

nodes (1,7). Similarly, nodes (1,7) are the interrelationship 
node pair for element corresponding pair ( (?) , © ); 
nodes (2,8) for ( (?) , © ) and nodes (3,9) for ( (?) , ©) 
and ( (7) , © ). Note that for any interrelationship 
node pair, one node is from the start substructure and the 
other from the repeat substructure. 
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The table below for Figure 2.9 summarizes the recommended steps 
indicated above. The interrelationship node pairs are shown 
encircled. 


Start 

Element Substructure 
Node 


Repeat 

Substructure 

Node 


End 

Substructure 

Node 




© 

© 

© 

© 



10 


11 

11 

11 


12 


As an additional example, see Figure 2.7(b). Following the 
steps above, the interrelationship node pair is found to be 
(2,6) for the zee stiffened panel. 
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2.6.4 Some Restrictions on Substructuring Non-Repetitive Panels 

Whenever a panel wi thout repetitive structural groups is 
idealized with more than one substructure type, care must 
be exercised in the substructuring process. The panel must 
be divided such that any one substructure contains only one 
unconnected piece of structure. To illustrate this, Figure 
2.10 is presented. This skin-stringer assembly is divided 
into three substructures. Node 3 is not attached to node 4, 
node 6 is not connected to node 7 and node 5 is used to attach 
the stringer to the skin. In Figure 2.10(a), the repeat and 
end substructures defy the restriction stated above, because 
in the repeat and end substructures, there are two discontinuous 
or unconnected pieces of structures. In the end substructure, 
the structure containing nodes 9 and 11 is not attached to the 
piece of structure that contains nodes 12 and 10, and in the 
repeat substructure, the piece of structure containing nodes 
7 and 8 is not attached to node 6. Therefore, there are two 
unconnected structural pieces in each of the two substructures, 
repeat and end. Figure 2.10(b) is an example of a valid 
selection of nodes for this problem. Each substructure contains 
only one continuous piece of structure. 

Another restriction is in the input order of the nodes. The 
nodes should be ordered such that the input list of nodes can 
be separated into distinct groups, each group related to a 
different substructure. The start substructure nodes should 
be input first, repeat substructure nodes, second, and the end 
substructure nodes last. What is important here is the input 
order of nodes, not the user node numbers. If the nodes of 
Figure 2.10(b) are input in the same order as the user node 
numbers, an incorrect data set will result. Nodes 9 through 
11 would be out of order. Instead, the input order should be 
(1 ,2,3,4,5,6,7,8,9,11 ,12,10). The first group of nodes for 
the start substructure would be (1 to 6), second group for the 
repeat substructure, (7,8,9,11) and the last group for the end 
substructure, (12 and 10). 
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Figure 2.10 
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2.7 


Upper and Lower Bound Loads 

The panel upper bound load ts required to increase the efficiency 
of the method used in extracting the lowest eigenvalue. This 
panel upper bound is def i ned as the least of the element upper 
bound loads. An element upper bound load is the buckling load 
of a plate element with interior sides completely restrained. 

For a plate element with specified boundary conditions on one 
side, the upper bound load is the buckling load for the plate 
with the interior side completely restrained and the exterior 
side subjected to the specified boundary conditions. 

The upper and lower bound loads can be specified as input to 
BUCLASP2. If the bounds are not specified, an approximate upper 
bound load will be computed (Reference 1) and the lower bound 
load will be set to zero. The minimum eigenvalue for the 
structure should be between these bounds. For most problems 
the upper bound load calculated as above works very well, but 
in some extremely rare instances it is possible that no eigen- 
values wi 1 1 be extracted within the bounds . For these cases , 
the calculated upper bound is automatically increased by 20 % 
and a search is made between the old and new upper bound loads. 

For a structure with buckling load very close to the upper bound 
load, it is possible that the minimum eigenvalue will not be 
extracted. If this occurs, the determinant values printed out 
for the guess loads should be inspected for any irregularities. 
Normally, this irregularity is evidenced by a steady decrease 
in determinant value followed by an increase as the load is 
incremented, with no indication of a root having been encountered. 
In such cases there is a very strong possibility for an eigen- 
value to exist between the load steps where the irregularity 
occurred. In such cases , these loads should be used as the 
bounds for the problem and re-run. 
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2.8 


Mode Shapes 


The relative displacements may be computed for the buckling 
load corresponding to the critical wave number. The relative 
displacements are computed for each element with respect to 
its local coordinate system. They can be computed at evenly 
spaced prescribed points along the width of any plate element. 
From the output of relative displacements, the mode shape can 
be hand plotted. 
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2.9 


Modelling of Riveted Connections 

This section deals with the modelling of riveted connections. 

As an example, see Figure 2,11. If the rivet is ignored in 
the idealization, the stiffener would be assumed to be bonded 
to the skin. A suggested modelling technique is presented 
here to more correctly depict the behavior of riveted connections. 

The rivet is idealized as a line at nodes 4 and 5. In Figure 
2.11(a) the plate elements (T) and (?) are offset at node 4 
as shown. Node 5 is equivalenced to node 4 to create an 
attachment between these nodes. The plate elements (F) , (?) , 

J) , and (?) are all joined at node 4, the rivet line, but 
node 2 of element(F)and node 3 of element(T)are not connected 
and they are free to deform independently. Nodes 6 and 7 of 
elements (?) and (?) respectively, are also free to move 
independent of each other. This modelling is approximate in 
the sense that the connection between the stiffener and skin 
is reduced to a line connection. 


First, make an initial run with this idealization shown in 
Figure 2.11(a) and investigate the mode shapes to see if they 
are feasible. For an example of mode shapes not being feasible, 
assume that the downward deflection of node 2 in the Z Q direction 
if greater than the downward deflection of node 3 in the same 
direction. Wherever this type of physical unreality occurs, 
join the elements together at these nodes and rerun the 
problem. For the example cited above, the connection of 
nodes 2 and 3 is accomplished by making node 2 equivalent to 
node 3. Offset element 3 as indicated in Figure 2.11(b). 
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In the original idealization shown in Figure 2.11(a) , nodes 
2 and 6 seem unnecessary for the first trial analysis. Plate 
elements (7) and (F) and similarly elements (7) and (F) could 
be combined into one element. The only reason for including 
nodes 2 and 6 in the example is to compare the displacements 
at these locations with nodes 3 and 7 of the stiffener, to 
see if any unrealistic deflections occur at these points. 
Actually, the displacements along the junction line between 
the skin and stiffener top flange should be investigated to 
check the feasibility of the deformed shapes. 
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(b) Reidealized Structure 


Figure 2.11 Riveted Connection Example 
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3.0 PROGRAM COMPUTER DETAILS 


3.1 Machine Requirements 

The BUCLASP2 program is written for the CDC 6600 series computers. It 
requires the use of a card reader, line printer, disk storage, a 
minimum of zero and a maximum of one tape drive. 

3.2 Operating System 

The program runs under the SCOPE 3.1 or KRONOS 2.0 operating systems. 
All system routines used are assumed to be standard CDC release. 

With the exception of four special purpose routine in COMPASS, all 
source routines are coded in CDC FORTRAN IV. The overlay loading 
feature is used. 


Two Langley Research Center routines are used. 

JPARAMS COMPUTER PROGRAMMING MANUAL 

Volume I, Section K5.1 date 6-9-71 

This routine makes certain job parameters available to the user. 

SIMEQ COMPUTER PROGRAMMING MANUAL 

Volume I, Section F4.1 date 8-1-68 

This routine solves a set of linear equations. 
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3.3 Storage Allocation 

The program will LOAD and EXECUTE with a field length of 77000g. 

The field length estimate for the program is dependent upon the 
data set. The space for the buckling determinant is dependent 
upon the size of the three blocks of the determinant definition. 

When the eigenvector and relative displacements are also wanted 
the user must take into account in field length allocation that the 
start substructure, each repeat substructure and the end substructure 
have to be stored in core. 

The data processing overlay prints an estimate of field length 
required for the other overlays. Table 3.1 gives field lengths 
for some of the test cases. 


3.4 Timing and Output Estimates 

Time consumption for one data set depends on various factors: 

a. Number of wave numbers investigated. 

b. Number of nodes, plate and beam element defined for panel. 

The size of buckling determinant is proportional to number 
of unconstrained nodes. 

c. Number of plate elements with stiffnesses which are the same. 

The actual computer time used for some of the test cases are 
presented in Table 3.2. 

In general, for single data set runs, and with the intermediate 
printout switches off, the program will aenerally generate from 
30 to 100 pages. The number of wave numbers is a primary factor on the 
amount of printout. The intermediate printout should be used with 
care as it can be voluminous. 


3.2 



3.5 Disk File and Tape Utilization 

BUCLASP2 uses a maximum of four files in addition to the standard 
input and output. The four files are scratch files used for internal 
information transfer. 


The files referenced are: 


(1) 

INPUT 

Standard Input (Card Reader) 

(2) 

OUTPUT 

Standard Output (Line Printer) 

(3) 

TAPE1 

Scratch File (Binary) 

(4) 

TAPE2 

Scratch File (Binary) 

(5) 

TAPE3 

Scratch File (Binary) 

(6) 

TAPE4 

Scratch File (Binary) 


3.6 Control Cards 

There are basically three modes of execution, from 

(a) source 

(b) relocatable binary 

(c) absolute binary 

In the following, use of specific control cards has been avoided, 
rather the required sequence to operations is specified. All file 
names with the exception of BUCLASP are arbitrary. All overlays have 
the name BUCLASP, thus a file BUCLASP is generated at load time. For 
the cases above 

(a) 1. Obtain a source file, PROG, from permanent storage 

(Cards, TAPE, permanent disk file, etc.). 

2. Compile source file placing relocatable binary on BPRG. 

3. Load BPRG. 

4. Execute from BUCLASP 

(b) 1. Obtain a relocatable binary file, BPRG, from permanent storage. 

2. Load BPRG. 

3. Execute from BUCLASP 
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(c) 1 . Obtain an absolute binary file, BUCLASP, from permanent storage. 

2. Execute from BUCLASP. 
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TABLE 3.1 


FIELD LENGTHS 


TA SET 

NUMBER 
OF BLOCKS 

HALF- 

BAND 

SIZE 

SIZE (IN 
OCTAL) FOR 
LOADING 

SIZE (IN 
OCTAL) FOR 
DISPL 

MAXIMUM NUMBER 
OF BLOCKS IN 
77 k 

2A 

7 

13 

60320 

33150 

68 

4A 

6 

9 

56610 

30700 

90 

5 

5 

17 

62040 

35350 

26 

6A 

7 

11 

60340 

33200 

50 

11 

6 

13 

61150 

35010 

27 

18A 

1 

48 

62620 

40510 



Assuminc 


Length of (1,0) Overlay 

(2.0) Overlay 

(3.0) Overlay 


45100 (Octal) 
56240 (Octal ) 
27350 (Octal) 
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TABLE 3.2 

COMPUTER TIME 



DATA SET 

ORDER OF 
BUCKLING DET. 

NUMBER OF WAVE 
NUMBERS 

TOTAL CP 
TIME 

2A 

60 

1 

14 

4A 

48 

1 

20 

5 

64 

1 

18 

6A 

72 

1 

17 

11 

100 

1 

26 

18A 

48 

1 

50 


Data Sets 2A - 11 were repeated substructure problems. 

Data Set 18A was a start substructure only. 

The CP time used varies nearly linearly with the number of wave numbers 

investigated 
DATA SET 

• 

NODES 

NUMBER OF 
PLATES 

BEAMS 

SUBSTRUCTURES 

2A 

7 

11 

0 

7 

4A 

8 

4 

3 

6 

5 

8 

10 

0 

5 

6A 

10 

9 

0 

7 

11 

16 

13 

2 

6 

18A 

12 

12 

0 

1 
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4.0 PROGRAMMED DIAGNOSTIC MESSAGES 

There are two types of error messages. The first kind is a message 
followed by routine name and number. This type is fatal to further 
processing* so an exit is taken. The second type is used in the 
input data processing routines to indicate a data error has been 
detected, but an attempt will be made to process the rest of the 
input data. This kind has a message, but no number. 


Routines which print diagnostic messages. 


DATAPRO 

FINDIN 

BLKDEF 

LCNTRL 

DM 

GAPUPP 

DB 

SBLKDT 

DBLERT 

STORE 

PRDMTX 

PROOT 

SOLVE 

DISPLAC 

FETCHD 

READTR 

FNLDIS 


Input Data Processing Routines 
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DATAPRO 


802 

TOO MANY 

NODES 


Card 4 

Field 1 

804 

TOO MANY 

ELEMENTS 


Card 4 

Fields 2,3 and 4 

806 

TOO MANY 

BEAMS 


Card 4 

Field 2 

808 

TOO MANY 

PLATES 


Card 4 

Fields 3 and 4 


810 BAD EQUIVALENT NODE DATA 

E Card 

812 BAD INTERRELATIONSHIP NODE DATA 

C Card 

901 NUMBER OF VALUES, III, IN M-LIST EXCEEDS 30 
Cards 6 and 7 

902 END OF FILE ENCOUNTERED BEFORE L 
L Card was not found 

903 NON ELEMENT DATA WAS ENCOUNTERED OUT OF ORDER 

904 NUMBER OF NODAL SPRING STIFFNESS SETS EXCEEDED NUMBER 
OF NODES IN DATA SET 

905 NUMBER OF PLATE ELEMENT OFFSET CARDS EXCEEDED THE 
NUMBER OF PLATE ELEMENTS IN DATA SET 
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906 NUMBER OF REPEAT SUBSTRUCTURE INTERRELATIONSHIP NODE 
PAIRS INPUT EXCEEDED MAXIMUM OF IIIII 

907 NUMBER OF EQUIVALENT NODE DEFINITION PAIRS INPUT 
EXCEEDED MAXIMUM OF IIIII 

908 INPUT MATERIAL PROPERTY TABLE FAILED 

909 INPUT OF THICKNESS TABLE FAILED 

910 NUMBER OF NODE DEFINITIONS INPUT, XXXXX, DID NOT AGREE 
WITH NUMBER SPECIFIED IIIII 

911 NUMBER OF ERRORS DETECTED IN PANEL SUBSTRUCTURE 
DEFINITION = IIIII 

912 NUMBER OF PLATE ELEMENTS INPUT EXCEEDED NUMRER SPECIFIED 

914 INVALID LOAD OPTION, IIIII, SPECIFIED 
Card 5 Field 1 

915 INVALID WAVE NUMBER SEARCH OPTION, IIIII, SPECIFIED 
Card 5 Field 2 

916 I II I ERROR(S) DETECTED IN DATA 

917 CURVED PANEL REPEAT SUBSTRUCTURE INTERIOR ANGLE EXCEEDS 
180 DEGREES 

Card 4 Field 6 

918 CURVED PANEL REPEAT SUBSTRUCTURE INTERIOR ANGLE MULTIPLIED 
BY NUMBER OF REPEAT SUBSTRUCTURES EXCEEDS 360 DEGREES 
Card 4 Field 6 
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919 CURVED PANEL REPEAT SUBSTRUCTURE INTERIOR ANGLE EQUALS ZERO 
Card 5 Field 6 

MATERIAL TABLE INPUT, BUT NO THICKNESS TABLE PRESENT 

THICKNESS TABLE INPUT, BUT NO MATERIAL TABLES PRESENT 

NUMBER OF MATERIAL TABLES INPUT EXCEEDS MAXIMUM OF III 

INVALID PLATE TYPE 
Card P Field 5 

NUMBER OF LAM I NAS EXCEEDS MAXIMUM OF 25 
Number of Plates Laminas <_ 25. 

Card P Field 7 

INVALID PLATE OFFSET SWITCH, III 
Card P Field 6 

INVALID MATERIAL PROPERTIES INPUT FORMAT OPTION, III 
Card P Field 8, or 

Card B Field 6 

INVALID N 22 LOAD SELECTION SWITCH, III 
Card P Field 9 

NUMBER OF RELATIVE DISPLACEMENT SUBDIVISIONS WAS NEGATIVE 
Card P Field 10 

INVALID MATERIAL TABLE INPUT OPTION, III 
Card P Field 11 
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INVALID BEAM TYPE 
Card B Field 4 

NUMBER OF LAYERS EXCEEDED MAXIMUM OF 35 
Number of Beam Layers £ 35 

NUMBER OF NODAL CARDS READ EXCEEDS NUMRER DECLARED 

ILLEGAL BOUNDARY CONDITION, IIIII ON NODE III 

NUMBER OF S CARDS, IIIII, DOES NOT AGREE WITH NUMBERS 
OF SPRUNG NODES IIIII 

USER ID, IIIII, ON S CARD, IIIII, INVALID 

USER ID, IIIII, IN REPEAT SUBSTRUCTURE INTERRELATIONSHIP 
PAIR, IIIII, IS INVALID 

USER ID, IIIII, IN EQUIVALENT NODE PAIR, IIIII, IS INVALID 
Card E 

USER ID, IIIII, IN, AAAAA, SUBSTRUCTURE DEFINITION IS INVALID 
Card D 

NODE! ( IIII) DOES NOT PRECEDE NODE2(IIII) OF, AAAAA, 
SUBSTRUCTURE DEFINITION 
Card D 

END NODE (IIII) OF, AAAAA, SUBSTRUCTURE IS NOT CORRECTLY 
RELATED TO START NODE(IIII) OF, AAAAA, SUBSTRUCTURE 
Card D 
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NODE(IIII) FOR PLATE(IIII) IS INVALID 
Card P Fields 3 or 4 


NODE ( 1 1 1 1 ) FOR BEAM (IIII) IS INVALID 
Card B Field 3 

USER ID, IIIII, FOR PLATE ELEMENT OFFSET SET, IIIII, 

IS INVALID 

Card F Field 2 

NUMBER OF PLATE ELEMENTS INPUT (IIIII) DOES NOT AGREE 
WITH NUMBER SPECIFIED 

NUMBER OF BEAM ELEMENTS INPUT (IIIII) DOES NOT AGREE 
WITH NUMBER SPECIFIED 

IIIIIIIIII, IS INVALID NODE NUMBER 

NODE I, IIIIIIIIII, IS NOT IN START SUBSTRUCTURE 

NODE J, IIIIIIIIII, IS NOT IN REPEAT SUBSTRUCTURE 

NODE, IIIII, IN START SUBSTRUCTURE SHOULD BE PART 
OF A REPEAT SUBSTRUCTURE INTERRELATIONSHIP PAIR 
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FINDIN 

FIRST NODE XXX GREATER THAN LAST NODE 

NO INTERIOR NODES BETWEEN FIRST NODE XXX AND LAST NODE XXX 
BEAMDF 

INVALID TABLE POINTER DATA 
See discussion on Table Input. 

PLTDEF 

INVALID TABLE POINTER DATA 
See Section 6.4. 
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BLKDEF 

708 I II II , ERRORS DETECTED IN SUBSTRUCTURE DEFINITION NODE 
PAIRS 

710 PANEL WITH REPEAT SUBSTRUCTURE BUT NO END SUBSTRUCTURE 
DEFINED 

711 SUBSTRUCTURE DEFINITION PAIR IS OUT OF ORDER 

713 NO SUBSTRUCTURE DEFINITION NODE INPUT 

714 PANEL MATRIX DEFINED IS NOT SQUARE 

715 START COLUMN DEFINITION OF SUBSTRUCTURE AAAAAA AND AAAAAA 
SUBSTRUCTURE IS NOT RELATED PROPERLY 

Examine Substructure Definition and Element Connectivity. 
See Section 2.6.4. 

716 END COLUMN DEFINITION OF SUBSTRUCTURE, AAAAAA, AND AAAAAA 
SUBSTRUCTURE IS NOT RELATED PROPERLY 

Examine Substructure Definition and Element Connectivity. 
See Section 2.6.4. 

717 ROW DEFINITION OF SUBSTRUCTURE, AAAAAA AND AAAAAA, 
SUBSTRUCTURE IS NOT RELATED PROPERLY 

Examine Substructure Definition and Element Connectivity. 
See Section 2.6.4. 

718 SUBSTRUCTURE AAAAAA AND SUBSTRUCTURE AAAAAA ARE NOT 
INTERRELATED PROPERLY 

Examine Substructure Defintion and Element Connectivity. 
See Section 2.6.4. 
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719 THE SUBSTRUCTURE DEFINITION OF START SUBSTRUCTURE AND 
ELEMENT CONNECTIVITY OF START SUBSTRUCTURE IS NOT CORRECT 
Examine Substructure Definition and Element Connectivity, 
See Section 2.6.4. 

720 NO EQUIVALENT NODES SPECIFIED FOR A PROBLEM WITH BEAM 
ELEMENTS 

721 IIIII, ERRORS FOUND IN PROCESSING EQUIVALENT NODES 

724 ELEMENT UNASSIGNABLE TO A SUBSTRUCTURE 

BOUNDARY CONDITION, II, NOT ALLOWED ON BEAM DEFINITION 
NODE, IIIII 

BEAM ELEMENT, IIIII WITH NODE, IIIII HAS NO EQUIVALENT 
NODE DEFINED 

EQUIVALENT NODES, IIIII, AND IIIII, HAVE CONFLICTING 
BOUNDARY CONDITIONS, IIIII, AND IIIII 


4.9 



LCNTRL 

901 DOUBLE ROOT COULD NOT BE SUCCESSFULLY BYPASSED 

902 ZERO DETERMINANT COULD NOT BE SUCCESSFULLY BYPASSED 

903 SEARCH FOR BUCKLING CONDITION DID NOT STABILIZE IN 
SPECIFIED NUMBER OF 'INTERAT.IO.NS 

Check for Panel Being in Buckled Condition for Zero Load 
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30 DOUBLE ROOT ENCOUNTERED IN UPPER BOUND CALCULATION 

6APUPP 

MATRIX REDUCTION FAILED 
DB 

906 MATRIX BLOCK DEFINITION FAULTY 

SBLKDT 

THE MATRIX HAS A ZERO ROW 

MATRIX SIZES OR RELATIVE POSITIONS ARE INCOMPATIBLE 
Zero Determinant at Block XXX 

DBLERT 

DOUBLE ROOT DETECTED 

STORE 

22 NODE NUMBER NOT FOUND IN SUBSTRUCTURE ROW DEFINITIONS 

24 NODE NUMBER NOT FOUND IN ‘SUBSTRUCTURE COLUMN DEFINITIONS 

PRDMTX 

5 ELEMENT NODAL INTERCONNECTIVE MATRIX SINGULAR 

PROOT 

ZARK FAILED TO CONVERGE IN THE MAXIMUM NUMBER OF ITERATIONS 
SPECIFIED 


ZARK FAILED - A ZERO IN THE PATH OF A SUBSEQUENT ONE 
Comnlex Root Found That is Not One of a Conjugate Pair 



SOLVE 

8 SINGULARITY ENCOUNTERED IN ATTEMPT TO PRODUCE STIFFNESS 

MATRIX 


DISPLAC 

INSUFFICIENT CORE IS AVAILABLE TO CONSTRUCT THE STIFFNESS 
MATRIX FOR THE EIGENVECTOR SOLUTION INCREASE FIELD 
LENGTH BY _____ 

THE STIFFNESS MATRIX CANNOT BE DECOMPOSED 


FETCHD 

14 ATTEMPT FETCH A NON-EXISTENT NODAL DISPLACEMENT COMPONENT 

READTR 

10 ID ON NODAL TRANFORMATION INVALID 

12 ID ON ELEMENTAL TRANSFORMATION INVALID 

FNLDIS 

20 ATTEMPT TO INVERT ELEMENTAL TRANSFORMATION MATRIX FAILED 
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5.0 Restrictions 

The following restrictions apply to this version of the program. 

5.1 Analysis Oriental Restrictions 

a. The special loading cases are described in Section 2.5. 

b. Each layer is orthotropic with respect to plate axes. 

c. Prebuckling deformations are ignored and linear theory (small 
deformations and linear elastic materials) is used. 

d. The cross-section is uniform in the axial direction 

e. Loaded edges of each element forming the cross-section are 
simply supported. 

f. The arc of a curved plate element has to be less than or equal 
to 180 degrees. 

g. Buckling load will not be computed for a structure already 
buckled by the input loads. 

h. When the buckling load is nearly equal to the upper bound load, 
it is possible that the buckling load will not be found. See 
Section 2.7 for a possible remedy for this problem. 

i. The curved panel repeat substructure interior angle (a©) must 
be less than or equal to 180 degrees. 

j. The total interior angle of a curved panel with repeated 
substructures must be less than or equal to 360 degrees. 
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5.2 Programming Oriented Restrictions 

a. Maximum number of elements definable in the three substructures 
is 35, and the maximum number of nodes is 35. 

b. Maximum number of types of plate elements is 10. These are of 
the same type in the sense that their lamina stiffness matrices, 
loadings and radii are of the same content. 

c. The maximum number of layers is 25 for plate elements. 

d. Maximum number of beam elements is 10. 

e. The maximum number of layers is 35 for beam elements. 

f. Two real roots of the determinant expression det (DT) = 0 are 
considered double if they differ by less than 0.00003%. 

g. The imaginary part of the complex roots of the determinant 
expression det (DT) = 0 is exactly set to zero if its numerical 
value is less than 10" 6 , or when it is less than 10~ 5 times the 
real part of the number. A similar test applies to the real 
part of the number. 

h. It is assumed that no coupling exists between bending and 
stretching when all of the elements of the B. . matrix are less 

■ J 

than 1.0. 

i. Maximum number of wave numbers that can be investigated in one 
data set is 30. 

j. The minimum number of elements is either 2 plates or 1 plate and 
1 beam. 

k. There must be at least 1 interior node in each block type. 

l. Node restrictions in blocks; see Section 2.10. 

m. An element cannot have nodes both in the start and end substructures, 
which means closed cylinder problems must be done with a single 
substructure. 
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5,3 Mathematical Oriented Restrictions 

In the buckling calculation the process is done strictly incore 
where the number of repeat blocks is arbitrary, but will be limited 
In practice by the accumulation of rounding errors. 

The eigenvector computation requires that the buckling load be well 
isolated. Problems which buckle without a determinant change of 
sign will have multiple eigenvectors for the buckling load. 
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6.0 PROGRAM INPUT 

The data input to this program consists only of cards, and no data tapes 
are required. 

6.1 Input Data Format 

Fields of 5 or less are right adjusted integers. Fields of 10 are 
floating point numbers with format FI 0.0. 

CARD 1 (8A10) 

Columns Column 1 must be blank. 

2-80 Title of run. Any characters anywhere in columns. This 

title is printed out in several strategic places in the 
output for the purpose of identification. 

CARD 2 (1615) Intermediate Print Control (I PC) Array 

The intent of the options on this card is for diagnostic checking. 

Columns 

I- 5 I PC ( 1 ) = 1 Print intermediate results (1) 

= Blank Suppress print 

This option prints out control lists obtained from input data. 

6-10 IPC (2) = 1 Print intermediate results (2) 

= Blank Suppress print 

This option prints out the P roots of the determinant 
expression det (DT) = 0 (equilibrium equations) during 
the upper bound calculation. 

II- 15 I PC ( 3 ) = 1 Print intermediate results (3) 

= Blank Suppress print 

This option prints out the elemental and merged stiffness 
matrices along with the determinant calculation information 
during the upper bound calculation. 


6.1 



Columns 


16-20 IPC(4) = 1 Print intermediate results (4) 

= Blank Suppress print 

This option prints out information about the element upper 
bound initial estimates and the actual upper bound. 

21-25 I PC ( 5 ) = 1 Print intermediate results (5) 

= Blank Suppress print 

This option prints out the p roots of the determinant 
expression det (DT) = 0 (equilibrium equations) during 
the panel buckling calculation. 

26-30 I PC ( 6) = 1 Print intermediate results (6) 

= Blank Suppress print 

This option prints out the elemental matrices during 
the panel buckling calculation. 

31-35 I PC ( 7 ) = 1 Print intermediate results (7) 

= Blank Suppress print 

This option prints out the merged stiffness matrices and 
determinant calculation information during the panel 
buckling calculation. 

36-40 I PC (8) = 1 Print intermediate results (8) 

= Blank Suppress print 

This option prints out the compacted stiffness matrix, 
the element p roots and the element w. . 
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CARD 3 

(1615) Program Control Array 

Columns 

1-5 

JPC(l) = 1 

Calculate panel upper bound only 


= 0 (Blank) 

Proceed as normal 

6-10 

J PC ( 2 ) = 0 (Blank) 

No relative displacements 


= 1 

Compute relative displacements 

11-15 

JPC(3) = 1 

Data check only is performed 


= 0 (Blank) 

Non data check only mode 

16-20 

JPC (4) 

Number of iterations allowed in 


'= Blank 

buckling calculation 

Number will be defaulted to 100 


= n 

where 1 <_ n <_ 100 

21-25 

JPC (5) 

Number of iterations used in eigenvector 


= Blank 

calculation 

Number will be defaulted to 4 


= T1 

where 4 <_ n <_ 10 


Option should be used with care. Can be used to find 
multiple eigenvectors for problems with coincident roots 
by running problem with n iterations and then with n + 1 
iterations. 
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CARD 4 Problem Characteristics 

Columns 

I- 5 Number of Nodes in the section (i.e. in start substructure* one 

repetitive substructure and the end substructure) 

6-10 Number of Beam elements in the section 

II- 15 Number of Flat Plate elements in the section 

16-20 Number of Curved Plate elements in the section 

21-25 Curved Panel Repeat Option 

= 0 (Blank) Non Curved Panel Option 

= 1 Curved Panel with multiple repeat blocks 

31-40 Curved Panel Repeat Block Included angle (a©) 

See Section 2.6.2 for discussion. 

41-50 Length of the section 

51-60 Transverse Load or Strain value 

61-70 Biaxial Load Ratio for first plate element 
Input only when this option is selected. 

For the transverse load or strain value, input a negative value for tension 
and a positive value for compression. Biaxial load ratios are computed with 
the same sign convention as above. 
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CARD 5 
Columns 

I- 5 

= 1 

= 2 
- 3 

6-10 MOPT 
= 1 

=2 

=3 

=4 

II- 15 MM I 

16-20 MMA 


Analysis Options 


Load Options 

N-j -j (force/length.) applied to first plate is varied with 
a N 22 transverse load of constant value applied to 
"affected" plates. The N 22 load can have an initial value 
plus be a linear function of as specified by biaxial 
ratio. (Generally* set N 22 = 0 for a biaxial ratio problem) 

N ;2 2 » transverse load is varied at the same value on all 
affected plates with an initial constant being applied 
to all plates. 

N 22 (force/length), transverse load is varied at the same 
value on all "affected" plates with an initial constant 
strain being applied to the first plate element. 

Option control for the loop on the longitudinal wave 
number M, 

Start the loop at wave number MMI and loop until a minimum 
load is found, then interrupt (Max. 30 loops) 

Start the loop at wave number MMI and loop to wave number 
MMA (Max. 30 loops) 

Start the loop at first value of M-list and loop through 
M-list until a minimum load is found, then interrupt 
(M-list must be input on cards 6,7) (Max. 30 loops) 

Start the loop at first value of M-list and loop through 
M-list. (M-list must be input on cards 6,7) (Max. 30 loops) 

Starting value for the loop on the longitudinal wave number 
M 

Set = 1 for options 3 and 4 

End value for the loop on the longitudinal wave number M 
Set = number of values in M-list for options 3 and 4 
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Columns 
21-25 NLW 

26-30 NUP 
31-35 

41-50 SLW 
51-60 SUP 

See Section 2 
discussion of 
criteria, see 


Lower Limit for root searching criteria 
(If left blank, NLW is defaulted to 0) 

Upper Limit for root searching criteria. If left blank, 

NUP is defaulted to 1 NLW » 0, NLW < NUP. The program 
does a complete solution of the data for each root searching 
criteria n/n+1 where NLW <_ n < NUP. (n is the 
number of roots below the trial load.) 

Default number of element subdivisions for relative 
displacement calculation (defaulted to 5, if left blank) 

Lower bound load (first plate element) for search interval 
(defaulted to 0.0, if left blank) 

Upper bound (first plate element) load for search interval. 
If left blank an upper bound is calculated for the structure 
for each M value. 

If specified SLW < SUP. Note that the same SLW, SUP 
pair is used for every M value to be searched. 

5 for discussion of load options and Section 2.7 for 
upper bounds. For a discussion of the root searching 
Section 6.1 of Reference 1. 
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CARD 6 

M- list (limited to 30 values) 
(List of wave numbers ) 

Columns 


1-5 

First value of M-l ist 

6-10 

• 

• 

Second value of M-l ist 

m 

75-80 

Sixteenth value of M-l ist 

CARD 7 

M-l ist (continued) 


(Use this card if necessary; if not omit) 

1-5 

• 

• 

Seventeenth value of M-l ist 

• 

65-70 

Thirtieth value of M-l ist 
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CARD T-l 


Lamina Thickness Table 


Columns 

I- 5 Letters "THICK" 

II- 80 ...» ...» value, .../ 

where value may be + nnn.nnn 
or +_ n.nnnn E + nn 

Number of significant digits is limited to 14 

The table is terminated by the character slash (/). If the entries of 
a table exceed one card, the table may be put on more cards using columns 
11-80 (columns 1-10 blank) with the last table value followed by a slash. 
The values in the table are delimited by commas and physical end of card. 

Omit card if Table option not used for element property input. See 
Section 6.4 for discussion of table use. 

Note: Only one thickness table is allowed per data set. 
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CARD T-2 


Lamina Material Properties Table 


Columns 

I- 5 Letters "TABLE" 

9-10 Table Number 

II- 80 value, .../ 

Same rules for this table as the one on CARD T-l 

Omit card if Table option is not used for element property input* 
See Section 6.4 for discussion of table use. 

Table Lengths 

Element Length 

Plate 4 

Rectangular Beam 3 

Circular Beam 2 
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CARD C 

Repeat Substructure Interrelationship Nodes 

Columns 

1-1 

Letter "C" 

3-5 

Node i 

Pair 1 (this card) 

8-10 

Node j 



Node 

i of start substructure is interrelated 

• 

• 

to Node .i of repeat substructure 

73-75 

Mode i 

Pair 8 (this card) 

78-80 

Node 3 



Reneat this card for each set of 8 interrelationship node pairs until all 
interrelationship node pairs have been input. 


Note : Interrelationship nodes must be interior nodes, and must be included 

if the number of reneat blocks >_ 2. 

See Section 2.6.3 for discussion. Omit this card when the number 
of repeat substructures is less than 2. 
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CARD E 


Equivalent Nodes 


Columns 

1-1 Letter "E" 


3-5 

8-10 


Node i 
Node j 


Pair 1 (this card) 


73-75 

79-80 


Node i 
Node j 


Pair 8 (this card) 


Repeat this card for each set of 8 equivalent node pairs until all 
equivalent node pairs have been input. See section 6.2 for 
discussion of equivalent nodes. 

Note : Equivalent nodes must be interior nodes. This card is used 

to specify those nodes which move together structurally. 
Omit this card if no equivalent nodes exist. 
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CARD N 
Columns 
1-1 
3-5 
10-10 


21-30 

31-40 


Nodal Definition 

Letter "N" 

User Node Number 
Nodal Boundary Condition 

= blank Interior 

1 Simple Support 

2 Clamped 

3 Free 

4 Sprung (Card S must follow) 

Y - coordinate of node in right-handed global coordinate 
system 

Z - coordinate of node in right-handed global coordinate 
system 

Internal node number is determined by input position, i.e. 
first node input is internal node one, etc. 
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CARD D Panel Substructure Definition 

Columns 

1-1 Letter "D" 

3-5 The number of substructures that the buckling determinant is 

divided into. This includes the start and end substructures 
plus the reneat substructures. 

8-10 User node number which the Lowest Internal Mode Number in 

Start Substructure 

13-15 User node number which is the Highest Internal Node Number in 

Start Substructure 

18-20 User node number which is the Lowest Internal Node Number in 

Repeat Substructure 

23-25 User node number which is the Highest Internal Node Number in 

Repeat Substructure 

28-30 User node number which is the Lowest Internal Node Number in 

End Substructure 

33-35 User node number which is the Highest Internal Node Number in 

End Substructure 

A substructure is left undefined by not specifying substructure definition 

node pair. 

Permissible Combinations 

Start Substructure Only 

Start and End Substructures Only 

Start, ReDeat (one or more) and End Substructures 


It is reconmended that the start substructure be larger of the 
three substructures for efficiency reasons in the determinant 
calculations. Note: See card N for definition of Internal 
Node Number. 
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CARD S 
Columns 


Nodal Spring Stiffness 


I- 1 Letter "S" 

3-5 User Node Number 

II- 20 Local w component for node 

21-30 Local 0 component for node 

31-40 Local v component for node 

41-50 Local u component for node 

See Section 2.4 for discussion of sprung nodes. 
Omit this card if there are no sprung nodes. 
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CARD F 
Columns 
1-1 
3-5 

11-20 0FF1 


21-30 0FF2 


31-40 0FF3 

41-50 0FF4 


Plate Element Offset 


Letter "F" 

User Plate Number 

Offset z q in the local z direction at the starting 
end (y = -b/2) of the plate element, measured positive 
in the positive z direction, from the negative 
surface of the element to the grid. 

Offset y in the local y direction at the starting 
0 

end (y = -b/2) of the element, measured positive in 
the positive y direction, from the end of the element 
to the grid. 

Offset z Q at the end y = +b/2 of the plate element. 
Measured similarly to OFF! . 

Offset v q at the end y = +b/2 of the plate element. 
Measured similarly to 0FF2. 

Omit this if no plate offsets exist. 

See Section 2.3 for discussion of offsets. 
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CARD P 
Columns 
1-1 


3-5 


8-10 

13-15 

20-20 


25-25 


29-30 

35-35 


40-40 


43-45 


50-50 


71-80 


Plate Definition 


Letter "P" 


User Plate Number 


Plate definition node I 
Plate definition node J 


| (user numbers) 


Plate element type 
= 1 Flat Plate 
2 Curved Plate 


Element Offset Switch 
= 0 (blank) No offsets for this plate 
1 Element has offsets 

Number of laminas in element 

Input format for element material properties 
= 0 (blank) P-1 or P-2 cards are used for materials 
property input 

1 Material properties are input through tables, 

N 22 load selection switch 

= 0 (blank) The N 22 load is applied this plate ("affected") 
1 This plate will have an effective N 22 load of ZERO 

Number of subdivisions used in relative displacement 
calculation. If left blank, the default number specified 
on Card 5 is used. 


Input option for material properties 
= 0 (blank) Engineering constants E,-| , E 22> v-| 2> G 
are specified 

1 Lamina stress-strain Matrix Q is specified 


12 


Curved Plate Element Radius (see Section 2.2.1) 
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CARD P-1 


Plate Thickness and Material Properties 

Columns 



1-10 

T 

Thickness of lamina 

11-20 

E 11 

E- modulus for direction 1 

21-30 

E 22 

E- modulus for direction 2. E 09 need not be entered 

L.C 

for isotropic laminas. 

31-40 

RNUA 

Poisson's ratio, v j2 

41-50 

^12 

G- modulus 


The 

subscript 1 denotes the longitudinal axis and the subscript 


2 the transverse axis of the plate local coordinate system. 
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CARD P-2 


Plate Lamina Thickness and Q. , Matrix 

* J 

Columns 

I- 10 Thickness of lamina 

II- 20 Q-ji element of lamina Q,. . matrix 

21-30 Q-j 2 element of lamina Q.j matrix 

31-40 Q 22 element of lamina Q. j matrix 

41-50 Q,,. element of lamina Q. . matrix 

^66 i j 

Q. . are the elements of the stress-strain matrix. 
U 

See equation (4.3) of Reference 1. 
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CARD B 


Beam Definition 


Columns 

1-1 Letter "B" 

3-5 User Beam Number 

8-10 Beam Definition Node (user) 

13-15 Type of Beam 

= 1 General 

2 Rectangular (may be laminated) 

3 Circular (may be laminated) 

18-20 Number of layers if laminated* <• 

25-25 Input format for element material properties 

= 0 (blank) B-l, B-2 or B-3 cards are used ror material 
property input 

1 Material properties are input through tables 

41-50 Angle between the local v-axis of beam element and the global 
Y-axis (measured clockwise, from the global Y-axis) 

51-60 Area of beam (Option 1) 

Beam element definition node must be equivalenced to a node 
of plate element. 
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CARD B-l 

Columns 

1-10 

11-20 

21-30 

31-40 

41-50 

51-50 

61-70 

71-80 


General Beam Properties 

E-modulus of beam element in longitudinal direction 

G-modulus of beam element material 

Moment of inertia about local v-axis 

Moment of inertia about local z-axis 

Warping constant for beam element about shear center 

Torsion constant of bead or lip about shear center 

v distance measured from the shear center to the centroid of 

the beam (measured parallel to local y axis) 

z distance measured from the shear center to the centroid 

of the beam (measured parallel to local z axis) 

See Section 2.2.2 for description of local coordinate system. 
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CARD B-2 


Rectangular Beam Lamina Thickness and Material Properties 


Columns 

I- 10 Lami na Thickness 

II- 20 E-modulus of lamina in longitudinal direction 

21-30 G-modulus of lamina material 

31-40 Width of beam element 

If the width of all lamina is the same, then only the width 
for the first lamina need be specified. 
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CARD B-3 Circular Beam Lamina Radius and Material Properties 

Columns 

I- 10 Lamina Outer Radius 

II- 20 E-modulus of lamina in longitudinal direction 

21-30 G-modulus of lamina material 

For a tube element, input the hole as the first lamina with zero E 
and 0. 
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CARD L 


Last Card of Data Set 


Columns 

1-1 Letter "L" 

This card must terminate each data set. 
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6.2 Use of Equivalent Nodes 

Equivalent nodes are used to cause two nodes to move together 
structurally. Beam are attached to plate elements so beam definition 
nodes must always be made equivalent to plate nodes. Care must be 
exercised so that (1) the use of equivalent nodes does not alter 
a substructure definition, and (2) a node that is a left (right) 
node of an equivalent pair does not appear as a right (left) 
node of a subsequent equivalent pair. See examples below: 


not permissable 


permissable 
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6.3 Data Stacking 


CARD 

1 

2 

3 

4 

5 

6 
7 
N 
D 

T-l 

T-2 


ORDER 

1 

2 

3 

4 

5 

Next if required 


C 
E 
S 
F 
P 

P-1 (or equivalent) 
or (not both) 

P-2 (or equivalent) 
B 


Omit any cards that are inapplicable 


Repeat for 
each lamina 


All plate definition cards 


B-l 

B-2 (or equivalent) Repeat for 
B-3 (or equivalent) each lamina 
L 


All beam definition cards 
(Omit if no beams) 

Last card in data set 


Repeat this order of card stacking for multiple data sets. 



6,4 Use of Tables 

The program has a table input option to facilitate the repetitive 
nature of the input for the thicknesses, radii and material properties 
of the plate and beam elements. This option replaces the P-1, P-2, 

B-2 or B-3 cards. The first item of the P-1, P-2, B-2 and B-3 cards 
(lamina thickness or lamina radius) is specified on the T-l card. 

All other items of the P-1, P-2, B-2 and B-3 cards are input on the 
T-2 card. These items are the material properties, E, v, 6 or Q, 
for each layer. Each set of values on the T-2 card will be identified 
by a user table number. In a data set, it is permissible to use both 
input options for different elements, but the same input option must 
be used for all layers of any one element. 

To use this table input, first input the proper parameter on columns 
35 and 25 of the P and B cards, respectively. Then, on subsequent 
cards, one or more for each element, specify the table numbers for 
element material properties and the locations of the thicknesses or 
radii on the T-l card. An example of the card is shown below. 

(1,2,3,,,,!/ 3,2, ,5, 6, 7,6) 

thickness material 

or radii properties 

The left and right parentheses are used to initiate and terminate the 
data for one element. The / is used to separate the T-l and T-2 
entries. The numbers to the left of / corresponds to the n^* 1 entries 
of the T-l card. The numbers are input in the order of the layers 
and separated with commas. If a number is omitted between commas, 
the previous integer is used for that layer. The numbers of commas 
to the left and right of / must always be equal. 
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The data for the example above will be interpreted in the following 
manner. 


Layer No . n 

1 

2 

3 

4 

5 

6 
7 


Item of T-l Card 

1 

2 

3 

3 

3 

3 

1 


Table No. of T-2 

3 

2 

2 

5 

6 
7 
6 
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7.0 OUTPUT 

The printing of the input data is coordinated with the description of 
the input in Section 6.1 . The items of output that are produced are 
as follows: 

1 . Header Page 

2. Program Control Options 
Card 3 

3. Problem Characteristics 
Card 4 

4. Analysis Options 
Cards 5, 6 and 7 

5. Nodal Data 
Card N 

6. Substructure Definition Data 
Card D 

7. Repeat Substructure Interrelationship Data 

Card C. Node I is in the start substructure and Node J is in 
the repeat substructure. 

8. Equivalent Node Data 

Card E. Node I is equivalent to Node J. 

9. Table Data 
Cards T1 and T2 

10. Element Data 

Cards P, P-1 (or equivalent) or P-2 (or equivalent). Each plate 
is accompanied by its A-, B- and D-matrices. The A-matrix represents 
the extensional stiffness of the plate element, while the B- and 
D-matrices are the coupling stiffness and bending stiffness, respec- 
tively. These matrices are computed with respect to the plate midplane. 

Cards B, B-l , B-2 (or equivalent) or B-3 (or equivalent) 

11. Plate Element Input Offsets 
Card F 

12. Plate Element Net Offsets 

The net Z0 offsets are measured from the midplane of the element 
to the point that is offset to. 
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13. Element Transformation Data 

The plate element widths, the sines and cosines of the transformation 
angles (local to global) at the definition nodes and the element 
types are printed. The sign convention for the angles are shown 
in Figure 5.2 of Reference 1. The plate width for curved elements 
is the arc length between nodes. 

14. Element Substructure Assignment Data 

Each element is assigned to a substructure by the following rules: 

a. If all nodes are in one substructure, then element is in the 
substructure. 

b. If one node is in the start substructure and one node is in 
the repeat substructure, then the element is in the repeat 
substructure. 

c. If one node is in the repeat substructure and one node is in 

the end substructure, then the element is in the end substructure. 

d. If one node is in the start substructure and one node is in 

the end substructure, then the element is in the end substructure. 
This assignment information is used in the total load calculation 
and the relative displacement calculation. 

15. Curved Panel Transformation Check 

The last node of the panel is transformed from its initial position 
to its final position using the number of repeat substructures and 
panel included angle. This is supplied for user visibility of the 
problem. 

16. Field Length Data 

The field length for the buckling calculation is exact, whereas the 
one for the eigenvector calculation is an upper bound to the one 
needed. 

17. Panel Upper Bound 

18. Buckling Load Search History 

The number of negative elements found on the diagonal of the panel 
matrix during the determinant calculation. The buckling determinant 

D 

- A* (2 ). There is a history printed for each wave number. The 
number of negative elements indicates the number of eigenvalues 
that are below the trial load. 
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19. Element Loads 

The inplane biaxial loads, and and the biaxial strains, 

E^-j and for plate elements are printed here. For beam elements, 
the axial loads and strains are printed. The total loads, axial 
and transverse are also included. Sign convention used here is 
such that tension is (-) and compression is (+). 

20. Buckling Load Computation Summary 

The critical wave number among those searched is identified. The 
buckling load for the first element of the panel is printed for 
each wave number. 

21. Eigenvector 

For each iteration, the eigenvectors for every freedom of all 
nodes, with the exception of constrained or sprung nodes, are 
printed along with the normalizing factors. User should inspect 
these eigenvectors to cheek on convergence. These values are 
printed with respect to the global coordinates. 

22. Relative Displacements 

The element displacements are relative to the element local 
coordinate system. 
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8.0 SAMPLE PROBLEMS 

The input, output and accompanying sketches of four (4) sample 
problems are included in this section. The circled items on the 
figures and discussions are user element numbers and the uncircled 
integers are user node numbers. 


8.1 Example 1: Potpourri Innut Problem 

The samole input for this problem is presented for this rather 
arbitrary structure, shown in Figure 8.1. This example is included 
here to illustrate many input features of BUCLASP2. 

The total structure is assumed to comprise of two more repeat 
substructures than is shown in Figure 8,1. Mode 101 is simnly 
supported, node 117 damned, node 107 sprung, and nodes 104 and 
113 are free to deform. Plates (Jt) and (J8) are separated by an 
offset of 1 inch. (16) is a three layered circular beam and (s) 
is a general beam element that is offset 2 inches from the midolanes 
of plate elements (6)and (jT) . 
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Figure 8. 1 Example No. 1 - Potpourri Input Problem 



POTPOURRI INPUT PROBUEH 
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17 

3 

13 2 0 



500. 

1 

2 

110 1 




n ioi 

i 

0.0 

0.0 



N 102 


1.0 

0.0 



N ICS 


1.0 

2.0 



H 104 


2.0 

2.0 



H 105 


5.0 

0.0 



H IQS 


15.0 

0.0 



H 107 

4 

21 .0 

-0.7 



N 106 


21 .0 

0.0 



N 109 


30.0 

0.0 



* 110 


36.0 

o.o 



W 111 


47.0 

0.0 



N 112 


60.0 

o.o 



N 113 

3 

60.0 

14.0 



N 114 


66.0 

0.0 



H 115 
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0.0 



N 116 
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2 
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0.0 
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F 7 
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.1 
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.3 

3.E7 

F 9 
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8,2 Example 2: Curved Sti ff ened Panel 


This example problem is illustrated in Figure 8.2, Rectangular beam 
elements (?) , (5) and (f|) represent local skin reinforcements and 
the integral stiffeners are idealized with flat plate elements (F) , 

(?) , (10) and (lj)l The sides of the panel are simply supported. 

N-ji critical load is sought for a constant transferse Toad of 
200 lbs. /in. 

This problem could be analyzed as a curved panel with repeated substructures. 
The use of this option could reduce this problem by several nodes and 
elements. 
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12 

3 



19.5549 

2.6477 


H 

13 

1 



21 .8305271 .2071442 


C 

4 

4 

7 

7 

10 10 



B 

3 

1 

6 

7 

10 11 

13 


F 

3 


0.0 


0.0 

.413 

0.0 

F 

5 


-413 


0.0 

0.0 

0.0 


7 


0.0 


0.0 

.413 

0.0 


9 


.413 


0.0 

0.0 

0.0 

F 

11 


0.0 


0.0 

.413 

0.0 


13 


.413 


0.0 

0.0 

0.0 


1 

1 

2 

2 

0 1 

-0 

-0 0 


0.163 


10.5 

E 6 


.32 

3.975 £6 


2 

2 

3 

1 

0 1 

1 

-0 0 


0.22 


10.5 

EE 


.32 

3.975 £6 


3 

2 

4 

2 

1 1 

~Q 

-0 0 


0.163 


10.5 

EE 


.32 

3.975 ES 


5 

4 

5 

2 

1 1 

-0 

-0 0 


0.163 


10.5 

EE 


.32 

3.975 £5 


6 

5 

6 

1 

0 1 

1 

-0 0 


0-22 


10.5 

E 6 


.32 

3.975 E 6 


7 

5 

7 

2 

1 1 

-0 

-0 0 


0.163 


10.5 

E 6 


.32 

3.975 £6 


9 

7 

6 

2 

1 1 

-0 

-0 0 


0.163 


10.5 

EE 


.32 

3.975 £6 


10 

6 

9 

1 

G 1 

1 

~G 0 


0.22 


10.5 

£6 


.32 

3.975 £6 


11 

6 

10 

2 

1 1 

-0 

-0 0 


0.163 


10.5 

E 6 


.32 

3.975 £6 


13 

10 

11 

2 

1 1 

-0 

-0 0 


0.163 


10.5 

EE 


.32 

3.975 E 6 


14 

11 

12 

1 

0 1 

1 

-0 0 


0.22 


10.5 

E 6 


.32 

3.975 £6 


15 

11 

13 

2 

0 1 

-0 

*0 C 


0.163 


10.5 

£6 


.32 

3.975 £6 


4 

4 

2 

1 



-3.62 


.3 


30. 

EE 

6 . £6 

.5 


B 

6 

T 

2 

1 



0.0 


.3 


30. 

E 6 

6. £6 

.5 


B 

12 

10 

2 

1 



3.62 


.3 


30. 

E 6 

6. £6 

.5 



L 


0 


196 * 

196* 

198. 

196. 

196. 

196. 

196. 

196 . 
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* * o t, % & n ton r a / B U CL ASP 2 cc u« 

8UCKUW* LOADS or LAMINATED STIFFENED PANELS 
LOADING — IN- PL A HE BIAXIAL 

lOJHBART CONDITIONS — CU 5IMPLT SUPPORTED ALONG EDGES PARALLEL TO THE TRANSVERSE AXIS 

<2J CLASTXCALLT RESTRAINED ALONG ANT EDGE PARALLEL TO THE LONGITUDINAL AXIS 
ttpes or elements — cn plat plate 

22) CURVED PLATE 

»> beam 

example no. 2 CURVED stiffened panel 


THIS IS AN EDITED EXAMPLE OF THE PROGRAM OUTPUT. 


8.8 



PROGRAM CONTROL OPTIONS 
BUCKLING LOAD AND RELATIVE DISPLACEMENTS 

NUMBER OF ITERATIONS ALLOWED IN BUCKLING CALCULATION = 100 

NUMBER OF ITERATIONS USED IN EIGENVECTOR CALCULATION S 4 

PROBLEM CHARACTERISTICS 


NUMBER OF NODES s 13 

NUMBER OF FLAT PLATE ELEMENTS = 4 

NUMBER OF CURVED PLATE ELEMENTS = 8 

NUMBER OF BEAM ELEMENTS = 3 

PANEL LENGTH = 50.000 

TRANSVERSE LOAD - N22 = 200.000 


ANALYSIS OPTIONS 


LOAD OPTION (LOADOP) = 1 

WAVE NUMBER SEARCH OPTION CHOFT) s 2 

INITIAL VALUE OF LONGITUDINAL WAVE NUMBER (MMI) = 1 

FINAL VALUE OF LONGITUDINAL WAVE NUMBER <MMA) = 3 

LOWER LIMIT FOR ROOT SEARCHING CRITERIA (NLW) = 0 

UPPER LIMIT FOR ROOT SEARCHING CRITERIA <NUP) = 1 
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« O 8 A L A r & 


COORDINATES BOUNDARY CONDITIONS 


SPRING STIFFNESSES 


NODE V 2 CODE MVO X (THETA) X(V) 


% 

-21.13053 

1.20714 

1 

z 

-1S.722Q0 

.$$710 

0 

3 

~t». 35480 

2.64770 

3 

4 

~12. 49170 

•39440 

0 

s 

-6.24900 

.09860 

0 

t 

-$.19310 

1.86630 

3 

T 

0.0QQQ0 

O.OGQOO 

0 

• 

$•24900 

.09860 

0 

» 

$•19310 

1.86630 

3 

10 

12.49170 

.39440 

0 

It 

IS. 72200 

.88710 

0 

12 

IS. 93480 

2.6477 Q 

3 

13 

21.83053 

1.20714 

1 

SUBS 

T R U C T U K 1 

i DEF1 Nl Tl 

C N 


WUWBE* OF SUBSTRUCTURES s 3 

SUBSTRUCTURE FIRST NODE LAST NODE 

start % $ 

REPEAT T SO 

end is 13 


C S U I VALENT NODE DATA 


NOCE I NODE J 

4 4 

T T 

so • 10 


R(U) 
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ccavfcp PLATE CLEMENT USES NUM2£E I USER NOC£ 1 S U$CR WOCC J 2 *A 5.2 Iff s 193.00000 

PLATE LOAD - N2 Z .* 2CQ .OCC 

NUMEE* or RELATIVE CISFLACC-MENT SUfi'P I VISIONS - 5 

material property input oft I CM s 0 

HATER I AL PROPERTY INPUT FORMAT OPTION = -0 


layer 

THICKNESS 

EXX 

EYY 

HUXT 

HUTX 

6 

1 

1 . STQQGOE-Ol 

1 . 059G^GE*37 

1 .Q<acaOE*0? 

9 . 2fl0G0GE-0i 

3.20QC0QE-01 

3.975C00E+06 


8.11 
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A-MATRIX 


190675 1 .337 
61 01 60 .428 
0.000 

B-MATRIX 

0.000 
0 . 000 
0.000 

c-matrix 

4221 .706 
1350.946 
0.000 


61 Gt 60. 428 
190675 i . 337 

o.ono 


o.ooo 


o.ooo 

0.000 


1350.946 
4221 .70S 
0.000 


0 

0 

647925 


0 

0 

Q 


0 

a 

143 V 


.000 

.000 

.000 


000 

000 

000 


coo 

OGQ 

560 
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PLAT PLATE ELEMENT USER NUMBER Z USER N^-CE ! Z USER HOSE J 3 

PLATE LCAC - HZZ z C.OCQ 

HUHacs or relative displacement *u?di v-.-rc^s = s 


HATER l At 

PROPERTY INPUT Of-riGH 


= 0 




HA TER X At 

PROPERTY INPUT FORMAT 

OPTION 

= -0 




LAYER 

THICKNESS 

zxx 

ETT 

Nuxr 

NUYX 

« 

1 

2 .2D0QQQE-Q1 

1..05GCCG£«07 

i • CSGC.3CE+Q7 

5.20CGGGE-0I 

■3.2.G5CC0.6-01 

3.973GQQZ+QS 
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A -MATRIX 


2573 529 .412 
823529 .412 
C. GOO 

B- MATRIX 

-.000 

-.000 

0.000 

C- MATRIX 

10379.302 
3321 .569 
0.000 


823529.412 

2573529.412 

0.000 


-.000 

-.000 

0.000 


3321.569 
1G379.902 
0 .000 


0.000 

O.OOQ 

874500.000 


0. O'jG 
0.000 
-.000 


0.000 

G.OGO 

3527,150 
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UAH tU.Mi.NT USCR NUMBER 4 USER NOCE 4 

MATERIAL property input format option 5 -o 

OPTION 2 

RECTANGULAR BEAM ELEMENT 


BEAM ELEMENT CONSIGNS 

WIBTH WB= .5000 

THICKfCSS TB= .3000 


BEAM ELEMENT PROPERTIES 

FOR LAMINATED BEAM ELEMENTS THE PROPERTIES AF£ GIVEN FOR ALL LAM I NAS 


AREA AFA = 
MOM. CF INERTIA RTY = 
MOM. OF INERTIA RZZ = 
POLAR MOM. OF ITCRTIA RIP - 
MAKING CONSTANT GAM = 
TORSION CONSTANT RJ = 
E-MOCULUS EXX s 
G -MODULUS GBA := 
Y DISTANCE YN = 
Z DISTANCE ZN = 


1 .5000E-01 
3.1250E-G3 
1.1250E-CB 
4.2500E-03 
2.3438E-G5 
2.8174E-03 
3.00Q0E+07 
8.0U00E+O6 
-8.88X8E-16 
0 . 
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F L A 

T £ i L £ M £ H T 

I H F U 

T 

OFF 

SETS 

EUEM. 

START 20 START TO 

ENC ZG 


END TO 


3 

Q.OOCOO 0.00000 

.41300 


0 .00000 


3 

.41300 0.00000 

0.00000 


0.00000 


1 

0.00000 O.OOGCO 

.41300 


G.OGGGO 


9 

.41300 O.OOGCO 

0.00000 


0.00000 


IV 

0.00000 0.00000 

.41300 


0.00000 


19 

.41300 O.OGCOQ 

O.OCGOG 


0.00000 


P L A 

T £ ELEMENT 

net 

O 

FFSE 

T 5 


ELEM. NO. 

START 30 

START TO 

END 30 

END 70 

1 

0.00000 

0.00000 

O.OGGCO 

O.OOGOO 

2 

0.00000 

0.00000 

O.OCGOG 

O.OGGGG 

3 

0.00000 

O.OOGCO 

.33150 

Q.OQCGG 

5 

.33150 

0.00000 

O.OGGGG 

O.OOGOO 

e 

0.00000 

0.00000 

O.OOGCO 

O.OGGGG 

7 

0.00000 

0.00000 

.33150 

O.CGGCO 

9 

.33150 

0.00000 

0.00000 

O.OGGGG 

10 

O.QGOOO 

0.00000 

O.OOGOO 

0.00000 

11 

0*00000 

0.00000 

.33150 

G.OGGGO 

13 

.33150 

0.00000 

O.OOGOO 

O.OGGGG 

14 

0.00000 

0.00000 

0.00000 

0.00000 

15 

0.00000 

0.00000 

O.OOGOO 

O.OGGGG 


8.16 



EtC 

M E 

NT T ft & 

n ft r c 

r ■* a r i c. 

N DAT 

k 


etc*. 

NO. 

WIDTH 

TRANSFOfit'ATlCN 



ct.TTPE 




SIMUT 

COS (I) 

SIN (4) 

CCS(rf) 



1 

5. 1249ft 

-.11 C26 

.9339G 

-.Gr»<56 

.93552 

CRVC FLA TC 


z 

l.fftftSft 

.99552 

.09434 

.99552 

.05434 

FLAT PLATE 


3 

ft. 25001 

-.09456 

.99552 

- .06309 

.99601 

CRV2 plate 


4 

ft. 2499 S 

-.06309 

. 998G1 

-.03156 

.99950 

CR VC PLATE 


3 

l.TftftSft 

.99950 

.03161 

.99950 

.C 3. 16.1 

FLAT PLATE 


ft 

ft. 25004 

-.C3155 

.99950 

.0GCG1 

l.OOCOO 

CRVC plate 


7 

ft. 2 5004 

--QCGQ1 

1 . GGOGU 

.03 105 

.99950 

CRVE PLATE 


ft 

1.7695S 

.99950 

-.c?:6i 

,99950 

-.03161 

FLAT PLATE 


ft 

ft. 2499ft 

.031 56 

.99950 

.06309 

.99601 

CR VC PLATE 


10 

ft .25001 

.06309 

.99801 

.09456 

. 99551! 

CRVC FL4TE 


u 

1.76852 

..99 i 5.2 

- . 09454 

.99552 

-.0345'. 

flat PLATE 


12 

3.1209 

.094 5*5 

.33552 

. X IC26 

.99390 

CRVC PLATE 


i 


0.C00C3 

1 ..CO-CC.O 



BEAK 


2 


C.OCOC*- 

i.OOGCO 



3EAM 


9 


.06214 

.99800 



BEAM 

E L E 

* E 

H T ft u e 

S T R U 

C T U R E 

A S S I 

S N M E N T 

DATA 


J*BC« 

Block 

1 

start 

2 

START 

3 

START 

4 

START 

5 

START 

ft 

REPEAT 

7 

S.CHSAT 

ft 

REPEAT 

9 

t 

10 

end 

11 

tNC 
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%Z END 

t START 

2 REPEAT 

3 REPEAT 
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at************************* 

ASSUMING A PROGRAM LENGTH OF 0000050000 OCTAL 

(EXCLUDING BLANK COMMON) 

FOR LOAD CALCULATION THE 

RECOMMENDED FIELD LENGTH IS 0000051430 OCTAL 

FOR EIGENVECTOR CALCULATION THE 

RECOMMENDED FIELD LENGTH IS 003CG52314 OCTAL 

************************************************** 
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MODE 


IS 


M = 


* * 3 * % * & $ $3 


EXAMPLE NO. 2 CURVES STIFFENED PANEL 


PANEL UPPER DOUND 


3.4521E+05 < ELEMENT 6 ) 


ABLE LOAD ON 

NUMBER C F NEG. 

BUCKLING DETERMINANT 

(A* 

T PLATE ELEfCNT 

ON DIAGONAL 

A 

8 

172602.9719 

9 

-.50675 

420 

86301 .48597 

8 

.24295 

412 

43150.74298 

7 

-.95983 

404 

21575.37149 

6 

.93331 

392 

10787.68575 

5 

-.48272 

392 

5393.842873 

4 

.12240 

392 

2696.921436 

3 

-.32735 

388 

1348.460718 

2 

.85340 

384 

674.2303591 

1 

-.13539 

384 

337.1151796 

0 

.09735 

388 

505 .6727693 

0 

.30105 

384 

589.9515642 

0 

.50746 

376 

632.0909617 

1 

-.08445 

384 

611.0212629 

1 

-.73550 

380 

600.4864136 

1 

-.37144 

380 

595.2189889 

1 

-.17482 

380 

592.5852766 

1 

-.07280 

380 

591 .2684204 

1 

-.33373 

376 

S90.C0999Z3 

0 

.08561 

376 

590.9392064 

1 

-.12437 

376 

590.7745993 

1 

-.31138 

372 

590.6922958 

0 

.52886 

372 

590.7334476 

0 

.10866 

372 

590.7540235 

1 

-.10138 

372 

590.7437355 

0 

.93021 

364 

590.7488795 

1 

-.78201 

368 

590.7463075 

1 

-.36194 

368 

590.7450215 

1 

-.15189 

368 

590.7443785 

1 

-.75020 

364 

590.7440570 

0 

.09010 

364 

590.7442170 

1 

-.33036 

364 

590.7441374 

1 

-.11967 

364 

590.7440972 

1 

-.24434 

360 

590.7440771 

0 

.60997 

360 

590.7440872 

0 

.17758 

360 

590.7440922 

1 

-.57048 

356 

590.7440897 

0 

.07321 

360 

590.7440909 

0 

.38524 

356 


COMKNTS 
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U t H L Ij 1 L *v A Z $ 


CUM. 


NO* 

NU 

N22 

EPSll 

EPS22 

P-BEAM 

% 

$90.74409 

200.00000 

3. 077 67 5 09E -04 

6 .40484423E-06 


Z 

710.94294 

0. 

3.07767509E-04 

-9. 84 85 6 027 E -05 


3 

590.74409 

200.00000 

3.07767509E-04 

6.40484423E-06 


4 

590.74409 

200.00000 

3.G77675G9E-Q4 

6.40484423E-06 


3 

710.94294 

0. 

3. 077675 09E -04 

-9. 84 85 6 027 E -05 


0 

590.74409 

200.00000 

3. 077 67 5 09E -04 

6.40484423E-06 


7 

590.74409 

200.00000 

3.07767509E-04 

6.40484423E-06 


S 

710.94294 

0. 

3.07767509E-04 

-9.84856027E-05 


9 

590.74409 

200.00000 

3. 077 67 5 09E -04 

6.40484423E-06 


10 

590.74409 

200.00000 

3 .077675 09E -04 

6.40484423E-06 


11 

710.94294 

0. 

3 .077675096-04 

-9. 84 85 6 027 E -05 


12 

590.74409 

200.00000 

3.07767509E-04 

6.40484423E-06 


13 



3 .07767509E-04 


1384.9538 

14 



3 .077675096-04 


1384.9538 

15 



3.077675Q9E-O4 


1384.9538 

total 

IOU 

AT BUCKLING 




AXIAL 

s 

3.5Q3E+04 




transverse 

a 

1 .00GE+04 





CLEMEN! 

TYPE 

CRVD PLATE 
PUT PLATT 
CRVD PLATE 
CfiVB iPLATE 
FLAT PLATE 
CRVD PLATE 
CRVD PLAT! 
FLAT PLATE 
CRVD FLAT! 
CRVD PLATE 
FLAT PLATE 
CRVD PLATE 
BEAM EL. 
BEAM EL. 
BEAM EL. 
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i-XAHPLC no.?, CURVEC STirFE!£: PANEL 


CUCKLINO LOU COMPUTATION $ U M M A R 


CRITICAL WAVE NUMBER M = 1 


M 

i 

Z 

3 


LOAD 

590.7440922 

1660.160301 

2465.235669 
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EXAMPLE H0,2 CURVED STIFFENED PANEL 


* EIGENVECTOR * 

$ AND * 

* RELATIVE DISPLACEMENTS * 

* * 

* M = | * 

* * 
Mc*******4***********4**4****** 


STIFFNESS MATRIX BANDWIDTH = 9 


eigenvector 


ITERATION NUMBER 1 

NORMALIZING FACTOR 6.3592S98E+G4 


SUBSTRUCTURE NODE 


W 


t*cta 


V 


start 


fCFEAT 


EW 


2 

4 

5 

7 

8 

to 

ii 


4 *48993 080E-01 
1 .OGOOCBOGE+GC 
6.52239149E-G1 
3.56400531E-G5 
-6.51228758E-01 
-9.986G9014E-G1 
“4.48G32545E-Q1 


1*3603 1333E-01 
1.45899131E-G2 
“9-70813322E-G2 
-t*06962690E-0l 
-9 .6933426 1EH32 
1 *45467582£“02 
1 .357777 56E-01 


9.861 12995 E -03 
5 .16335990E-02 
4 .30934 805E-G2 
6 . 89987 954E-02 
4.306820G8E-Q2 
5 .16106616E-G2 
9. 86467 5 26 E “03 


EIGENVECTOR 


ITERATION NUMBER 2 

NORMALIZING FACTOR 3 .4.22 9813 E +07 


SUBSTRUCTURE NOCE 


W 


THETA 


V 


START 


JCFEAT 


END 


Z 

4 

5 
f 
8 

10 

11 


“4 .48796704E-CU 
-l.OOOOOOOOEHa 
-6.52157OS7E-01 
1 *3957085 9E “04 
6.52229265E-01 
9* 9977 6 003 E -Cl 
4 .48649627E-01 


-1 *36007 908E-01 
-1 .45593499E-02 
9 - 70601 762E “02 
1 .07097991E-01 
9. 7006034 9E -02 
-1 *45867 H9E -02 
-1.35966246E-01 


-9.87724972E-03 
-5 *167 346 15 E “02 
-4.31174909E-Q2 
-6 .90690919E-C2 
*4 . 3 1 12965 3 E “02 
-6 . 16461 913E-02 
-9. 866 5 896 IE “03 


EIGENVECTOR 


ITERATION NUMBER 3 

NORMALIZING FACTOR 3 .4250420E+07 


SUBSTRUCTURE NOCE 


TTCTA 


V 


START 


REFEAT 


2 

4 

5 

7 


4*48796704E-01 
1 .OOOOOOOOC+OG 
G.52157087E-01 
-1*3957 086 OE -04 


1 .36007 908E “01 
1.4 5 5 934 99E -02 
-9.70601762E-02 
“1.07097991E-01 


9.87724972E-03 

5.16734615E-02 

4.31174909E-O2 

6.90690919E-Q2 


U 

6 *204 832 18£ -03 
“8 *34 83604 6 E “03 
1-04 17027 9E -02 
2*5677 S098E-Q5 
-1-G3561897E-Q2 
8 -374235 85 E “03 
“6-170233 94E-03 


U 

“6-132523896-03 
8*372201 5 2E “03 
“1 - 03 95 3037 E -02 
“1 -66525 093E-06 
1-C3 96366 9E -02 
-8-36 885 93 8E -03 
6.1 81284 60E -03 


U 

6 *162523896-03 
“8 -37 2201 5 2E “03 
1.C3953G37E-02 
1 -66525G95E-G6 


8.23 



6 -6.522292652-01 -9 .7 006 034 3E -02 4,311296332-02 -1 .033636692-02 

10 -9.937760032-01 1.458671192-02 5.1646 1 9132-02 8.36 66:3382-03 

CHE II -4 .486496272-01 1.359662462-01 9. v6658961E-03 -6.181284 802-03 

EIGENVECTOR 


ITERATION number 


4 




NORMALRING r actor 


3.42504202+07 




SUBSTRUCTURE 

NQCE 

W 

THETA 

V 

U 

START 

2 

-4.487967042-01 

-1 .36GG7908E-01 

-9.87724972E-03 

-6 •18252389E-03 


4 

-1 .000000002+00 

-1 .455934 99E-02 

-5. 167346 15E -02 

8.37 2201 52E "03 


s 

-6.521570872-01 

9.T0601T62E-02 

-4 .31 174 909E -02 

-1 .039530372-02 

REPEAT 

T 

1.395708602-04 

1 .07097991E-01 

-6 .906909 19E-02 

-1.66525095E-06 


e 

6.522292652-01 

9. 7 006 034 9E -02 

-4 .31 1296 5 3E -02 

1. 0396366 9E -02 


10 

9.997760032-01 

-1.45867119E-02 

-5.16461913E-02 

-8. 36 885 93 8E -03 

ESC 

u 

4 .486496272-01 

-1 .35966246E-01 

-9. 866 5 896 IE -03 

6.1 81284 80E -03 
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EXAMPLE m«Z CURVED STIFFENED PANEL 

# ® 

* RELATIVE DISPLACEMENTS * 

# ^ 

* m 

* W 2 1 * 

* * 
^***^*$*******j^«jM£4******** 

START SECTION —BLOCK NO, 1 


ELEMENT NO.= 1 TYPE = CRVD PLATE NOCE 


Y 

W 


V 

-1.5625 

-.00000 


•03636 

-.7812 

-.11461 


.03612 

.0000 

-.22813 


•03541 

.7812 

-.33948 


.03424 

1 .5625 

-.44772 


.03261 

EL£*CNT NO . 2 

2 TYPE 2 

FLAT PLATE NOCE 

Y 

W 


V 

-.8843 

-.03260 


-.44772 

-4421 

-.09263 


-.44772 

0.0000 

-.15250 


-.44761 

-4421 

-.21230 


-.44739 

.8843 

-.27208 


-.44706 

ELEHlNT NO. = 

3 TYPE = 

CRVD 

PLATE NOCE 

Y 

w 


V 

-5.1250 

-.44772 


.03260 

-1 .5625 

-.65203 


.02806 

.0000 

-.82402 


.02198 

1.5625 

-.94453 


.01471 

3.1250 

-1.00127 


.00670 

ELEICNT NO. 2 

5 TYPE 2 

CRVD 

PLATE NOCE 

Y 

W 


V 

-3.1250 

-1.00127 


.00669 

-1.5625 

-.99079 


-.00151 

0.0000 

-.91896 


-.00939 

1.5625 

-.79959 


-.01649 

3.1250 

-.65319 


-.02251 

ELE&CNT NO . 2 

6 TYPE 2 

FLAT PLATE NOCE 

Y 

W 


V 

-.8843 

.02248 


-.65319 

-.4421 

.06532 


-.65320 

0.0000 

.10606 


-.65305 

.4421 

.15073 


-.65275 

.8843 

.19339 


-.65228 


ELEMENT NOoS 4 TYPE s BEAM EL, NODE 


I = l NOCE J = 2 WIDTH 2 3.1250 

nCTA 

-.14676 

-.14605 

-.14396 

-.14056 


I = 2 NODE J = 3 WIDTH 2 1 .7685 

T>CTA 

-.13601 

-.13557 

-.13530 

-.13520 

-.13526 

I = 2 NOCE J 2 4 WIDTH = 6.2500 

T*CTA 

-.13601 

-.12262 

-.09521 

-.05753 

-.01456 

I 2 4 NOCE J = 5 WIDTH 2 6.2500 

THETA 

-.01456 

.02739 

.06301 

.08744 

.09706 

I r 5 NOCE J 2 6 WIDTH 2 1.7686 

T>€TA 

.09706 

.09675 

.09656 

,09648 

.09653 

= 4 
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w 

i B 00000 


V 

05167 


THETA 

01456 


REPETITIVE SECTIOH —BLOCK HO. 2 


ELEMENT N©.v ? TYPE 5 CRVD PLATE NODE 


Y 

W 

V 

-3.1250 

-.65319 

-.02252 

-1.5625 

-.49764 

-.02731 

0.0000 

-.33523 

-.03078 

1.5625 

-.16816 

-.03267 

3.1250 

•00014 

-.03357 

ELE»CNT NO. = 

9 TYPE s 

CRVD PLATE NODE 

Y 

W 

V 

-3.1250 

.00014 

-.03357 

-1.5625 

.16843 

-.03287 

o.ocoo 

•33545 

-.03077 

1.5625 

.49779 

-.02730 

3.1250 

.65326 

-.02251 

ELOCNT NO.- 

10 TYPE = 

FLAT PLATE NODE 

r 

W 

V 

-.8843 

.02248 

.65327 

—*4421 

.06529 

.65327 

0.0000 

.10800 

.65313 

*4421 

.15065 

.65282 

.8843 

.19329 

.65235 

ELEJCNT NO.- 

11 TYPE = 

CRVB PLATE NODE 

Y 

W 

V 

-3.1250 

.65327 

-.02251 

-1.5625 

.79958 

-.01648 

0.0000 

.91886 

-.00938 

1.5625 

.99062 

-.00151 

3.1250 

1.00104 

.00670 

ELEJCNT N0.= 

8 TYPE = 

BEAM EL. HOCE 


W 

V 


.00014 

-.06907 

ELEICMT N0.= 

12 TYPE = 

BEAM EL. NODE 


W 

V 


1.00104 

.01158 

EMC SECTION —BLOCK HO. 3 

ELEICNT M0.= 

13 TYPE = 

CRVD PLATE NODE 

r 

w 

V 

-3.1250 

1.00104 

.00670 

-1 .5625 

.94427 

.01471 

•0000 

.62377 

.02198 

1.5625 

.65182 

•02806 

3.1250 

•44757 

.03260 


I = 5 NODE 4 s ? WIDTH = 6,2500 

THETA 
.09706 
.10175 
•10559 
•10756 
• 10710 

1= 7 NODE 4 s 8 WICTH s 6-2500 

THETA 

.10710 

.10754 

.10556 

.10171 

.09701 

I = 8 NODE J = 9 WIDTH = 1-7686 

T*€TA 

.09701 

.09670 

.09650 

.09643 

.09647 

I = 8 NODE J = 10 WIDTH = 6.2500 

theta 

•09701 

.08739 

•06295 

.02735 

-.01459 

= 7 

TVETA 

.10710 

= 10 

THETA 

-.01459 


I = 10 NODE J s a WIDTH = 6-2500 

TfrCTA 

-.01459 

-.05754 

-.09519 

-.12259 

-.13597 
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CLTMrnt tiO.t 14 ITA*. S rtAT n&TE WCt£ 1 s It HOtt 4 = It WltlH : 1-^6® 


¥ 

w 

¥ 

THETA 

• . §943 

.01012 

• t4m 

-.04000 

-•4421 

* . 00TS4 

.147*3 

- .03907 

0.0000 

* .0241S 

.14710 

-.03979 

. 44*1 

- .0421] 

.14700 

-.03973 

-••43 

* .0(030 

.1409* 

-.03970 


CLCWtMT NO.* 13 TT7E s C*V® OCATC WOOC I * 11 NODE 4 = 13 WIDTH = 3.1230 


▼ W ¥ THETA 


-1.90*9 

.14720 

-.01000 

-.04000 

-.▼•1* 

.11421 

-.01029 

-.04459 

.0000 

.07000 

-.01307 

-.04007 

•7012 

.03937 

-.01502 

-.09024 

1.3029 

-.00000 

-.01533 

-.05090 
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timing summakt 


FUNCTION CP TIME PER CENT OF TOTAL 


data processing 

.77 

.009 

UPFER BOUND CALCULATION 

14.04 

.156 

STIPFNE3S GENERATION 

71.90 

.608 

determinant CALCULATION 

1.61 

.018 

ELEMENT DISPLACEfCNTS 

.63 

.007 
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8.3 Example 3: Advanced Structural Panel 

The geometry and idealization for the panel is shown in Figure 8.3. 
The user should pay particular attention to the signs of the 
radii of curvatures of the curved plate elements. 
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Geometry 

Material Properties 

E (lbs/in 2 ) 10.5x10® 

v 0.32 

G (lbs/in 2 ) 3.975x10® 


Figure 8.3 Example No. 3 • Advanced Structural Panel 
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t*A*HE NO. 3 ADVANCED STRUCTURAL PANCt, 
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8.4 Example 4: Flat Stiffened Panel 

Figure 8.4 contains the pertinent data for an angle stiffened, bonded 
titanium panel with local boron reinforcements on the stiffeners. 

The Yg axis passes through the midplane of element (T), which is 
idealized as three bonded layers, skin, adhesive and stiffener flange. 
The boron reinforcements, beam elements (T) and (T), are idealized 
with 15 layers each. The boron reinforcement is actually simulated 
with six layers and the layers are separated by non-structural elements 
to retain the correct lamina geometry including the adhesive thicknesses 
This representation should give the proper stiffnesses for the boron 
reinforcement. 

Use is made of the flat panel repeat substructure option to minimize 
the input and also the computer costs. 
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Idealization 


Material Properties 



Figure 8.4 Example No. 4 — Flat Stiffened Panel 
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9.0 PROGRAM DESCRIPTION 


This section is a description of the organization and function of the 
various routines included in BUCLASP2. 


9.1 Overlay Structure 

BUCLASP2 consists of a (0,0) level overlay which acts as a monitor in 
selection of primary overlays. Each primary overlay performs a discrete 
logical step in the analysis. The overlay structure is: 


(1 ) (0,0) Overlay 
(a) BUCLASP, 0,0 

(2) Primary Overlays 

(a) BUCLASP,! ,0 

(b) BUCLASP, 2,0 

(c) BUCLASP, 3,0 


S0337A, Monitor 


DATAPR0, Reads all user input and outputs 
basic job description 

LOADING, Generates element stiffness, merges 
panel stiffness matrix, controls calculation 
of buckling load and calculates determinant 
of panel stiffness matrix 
DISPLAC, Calculates for buckling load panel 
eigenvector and element relative displacements 
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Main Overlay 0.0 

Program Name: S0337A 

Subroutines: ERROR 

PRINT 
I GAL 
PAC 
UNPAC 


Primary Overlay 1.0 


Main Prog: DATAPRO 

Subroutines: RDVLS 
BEAMDF 
USERID 
BLKDEF 
PLTDEF 
TBPOINT 
NEXTC 
FIND IN 
CL NODE 
RDTBLE 
IBLKNO 
STRMOV 


Primary Overlay 2.0 


Main Prog: LOADING 

Subroutines: LCNTRL 
DM 

UPPRBD 

GAPUPP 

TRANF 

SBLKDT 

BLKGEN 

SELIM 

REDUVE 

DBLERT 

STORE 

PRDMTX 

PROOT 

FPLATE 

FPLCOL 

BEAM 

SOLVE 

DT 

RGEN 

CDTM 

VIPDR 

ZARK 

MATZ 

TRANS F 

EVALUE 


Primary Overlay 3.0 


Main Prog: DISPLAC 

Subroutines: BANDW 
COMPAC 
BLKGEN 
SYMBND 
REDUVE 
DIS 
FETCHD 
READTR 
FNLDIS 
TRANW 
TRAN1 
FANGLE 
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LOADING 


Calculate 
Upper 
Bound 
\ ? ^ 


UPPRBD 

Calculate Panel 
Upper Bound 


LCNTRL/PB 

Form Stiffness 
Matrix and 
Determine 
Buckling Condition 


^ AH 'N 
M's Been 
Done 
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9.3 Program/Subroutine Description 
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13. L cards F-22 

14. Beam element F-23 

15. Flat Plate element F-23 

16. Curved Plate element F— 21 

User Diagnostics F-40 - F-90 
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b) Start and End substructures only F-38 

c) Start, Repeat and End substructures F-22 

Problem with 

a) Plate with zero B-matrix F-29 

b) Plate with non-zero B-matrix F-7 
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EVALUE The routine EVALUE extracts F-17 

the lowest eigenvalue from a 
general eigen problem of up 
to order 15 for the GALERKIN 
method. 



FANGLE For an element the local to 1. Non-curved panel option 
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For the last node in a sub- 
structure definition pair 
find the largest interior 
internal node smaller than 
the last node. 
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Zero determinant encountered at 
trial load 



LOADING LOADING controls the computation 1. Various wave number search options 
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remaining bits of A are unchanged. 
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up to 8 columns at a time and all 
rows. 



CD 

E 


1 — 
oo 


r>- 

CM 

CM 

CO 

T3 

P 

=3 

£ 

JD 

in 

P 

to 

r— 

LU 

i 

i 

s 

i 

E 

3 

QJ 

1 

| — . 

U- 

Lu 

LU 

Lu 

Id 

to 

P 

LU 



X 
















«r— 










O 






S- 










O 






P 





P 


P 



p 






id 

X 




E 


e 









E 

•I*" 

-a 

-O 


rd 


rd 



* 







L. 

CD 

O) 


U 


O 



s- 






CO 

P 

P 

to 


•r— 


•r** 


jz 

CD 




to 



rd 

rd 

3 


P 


P 


to 

p 

-a 





o 

E 

CD 



•r* 


*r 


rd 

U 

CD 





s- 

1 

CL 

to 


E 


E 


i— - 

rd 

“O 



>> 


0) 

CO 

CD 

"O 


cn 


cn 


to 

i- 

CD 



rd 


N 


L. 

u 


•r* 


•r- 



rd 

CD 



U 


| 

o 


rd 


to 


to 


cn 

x= 

O 


-a 

u 


e 


to 

o 






e 

o 

X 


CD 

rd 


o 

CD 

•1— 



<=r 




•r— 


<D 


u 



e 

N 


■i— 


r— 


i— 


to 

TJ 



<d 

cn 




-c 

<d 






to 

rd 

P 


p 

E 


x: 

-E 

u 

u 


Jr 


JC 


*! — 

jo to 

•f— i 


E 



p 

-P 

•i— 

•i— 


P 


p 


E 


E 


3 

> 


•»— 

•r* 

JZ 

to 


*p~ 


*r— 



« a> 

■r— 

p 

O 

•r— 


5 

5 

£ 



3C 


5 


•» 

rd P 

i — 

3 

o 

CD 





-C 






CD 

P o 


CL 

E 

a 


P 

P 

-P 

CL 

rd 

rd 

rd 

rd 



rd rd 

CD 

e 

CD 

CD 

OO 

E 

E 

E 


P 

P 

P 

P 


■r- 

“O t- 

N 

*r- 


S- 

O 

cy 

CD 

<u 

CD 

<d 

rd 

rd 

rd 

to 

P 

rd 

*r~ 


CD 


i — i 

E 

E 

£ 

S- 

-a 

"O 

"O 

-a 

P 

I 

r-.JC 

to 

rd 

r— 

P 

1 — 

<u 

CD 

CD 

o 





*r— 

P 

ro a 


P 


o 

CO 

r— 

r— 

p« 

E 

p 

P 

p 

p 

X 

O 

cn 

<D 

rd 

P 

x> 

1— 1 

CD 

CD 

CD 


rd 

rd 

rd 

rd 

CD 

1 

CD >> 


■a 


E CD 

DC 




L- 

E 

E to 

E 

E to 


~0 

r- E 

jQ 


P 

3 T3 

LU 

CD 

CD 

CD 

o 

S- 

P 

SL 

S- P 

L. 

.E 

i— rd 

id 

o 

o 

E CD 

h* 

P 

P 

-P 


,o 

O -r- 

,o 

O T- 

O 

CD 

*r- E 

p 

E 


•r- CD 

O 

1 0 

rd 

rd 

o 

P 

P a> 

P 

p cn 

z~ 





-o 

X u 

C 

i— 

r-~* 

r — 



■r— 


•r~ 

L- 


,•*— s. 

„ — 

, 

E 

rd X 

OC 

CL 

OL 

CL 

H- 

LU 

LU “O 

lu 

LU “O 

LU 

rd 

-Q 

a 

-a 

LU 

^ CD 

c 
















=c 

• 

• 

• 

• 

• 

* 

» 

• 

• 





• 

• 

o 

« — 

CM 

oo 

i— 

CM 

CO 

*3* 

LO 

to 






CM 



P U 

to ~o 

EE P 

E rd <d O 

3 

r~ C tO XJ P 


o p 

O rd » — CD •»— • 


P 

u g CD x: E to 


to 

- ' rd E •i-'* P *r* P 


P to 

P O P* ! e r— 


O E 

to rd o to cn 


O O 

CD TD W *r • *f- 


-r- 

cn r d -o 


P 

rd P P CD 


(D rd 

E c -a p p 


«C 3 

-|“ £ S- ^ CD E 


p cr 

o ro E <d* 


a» 

XJ CLP U JE -r- o 


to • 

i- Id tO r™* •r— 


CD E P 

Id o> E r— Id OJ P 


P 3 E 

U C C Id r— x> ‘ 1 - 


rd *r- <D 

•i- o a to e 


»— S- E 

E p p *r- TJ as 


3 JD CD 

O Id to • *r- 


U -r- r— 

p o “O >> to ro to 


r— r— <D 

Pr-r-r £~ O 


rd •*— 

P CD CL CD E 

LU 

O 3 CD 

cn •»— p r-^ a> 

CO 

CTP 

E ^.p P CD fd CD 

O 

P CD rd 

•1-0 1 O E o P 

Ol. 

o «— 

“O 00 CD T- L. 

DC 

O CD CL 

id i a) *o f— o> 3 

ZD 

DC SZ 

CD r— S- E CD O O 

CU 

CU P rd 

QC r— P OJ*Or-P 


LU 


O 

cc 


o 

§ 

D- 


CQ 


Q 

C£ 


CO 


Q 

q: 


9.19 



) — 







cd 

00 

*— 


o 


CM 


U1 

r— 

f— 1— - 




CO 

CM 

CM 

CO 


ro 


ro 

ro 

LU 

1 

1 1 


i 


i 

§ 

i 

1 


8 


i 

i 

f— 

LU 

Lu Lu 


Lu 


LU 

LU 

Lu 

LU 


LU 


Lu 

LU 




X 






X 









•r- 






•r— 









S- 






u. 









+3 















r 0 





X 

rd 





X 



"O 

E 


X 



•r* 

E 

X 




*r- 



0) 

i 


•i— 



S- 

1 

•r* 




u. 



+3 

CQ 





+3 

CQ 





-M 



</> 



+3 



rO 


+3 




rd 


"O 

QJ 

O 


rd 


X 

E 

O 

<d 


X 


E 


0) 

3 

U 


E 


•r— 

l 


E 


«r— 


l 


+3 

cr 

<u 


i 


U. 

CQ 

QJ 

i 


S- 


CQ 


to 

QJ 

N 


CQ 


+-> 


N 

CQ 


+J 




qj 

C- 





ro 

O 




rd 


O 


3 


-C 


O 


E 

L- 

SZ 

O 


E 


c- 


cr 

c 

4-> 


C- 


1 

0) 

+3 

u. 


l 


QJ 


0) 

o 

•r- 


CU 


CQ 

N 

•i— 

QJ 


CQ 


rsi 


c. 

•r* 

* 


N 



1 

* 

N 




l 



+3 



1 


O 

C 


1 


O 


c 


c 

rd 

+3 


c 


i- 

o 

+3 

cr 


u 


o 


o 

g 

CT 


o 


QJ 

c 

sr 

o 


QJ 


c 



l! 

QJ 


c 


N 


QJ 

c: 


N 




+3 

o 

E 





jz 

E 





-C 


«3 

4— 

0 ) 


JZ 


JC 

+-> 

QJ 

-C 


JZ 


4-> 


£ 

to 

r— 


+3 


+3 

•r- 


+3 


-M 


•i— 



cr 

LU 


•r— 


*r— 

5 

LU 

■r- 


•i— 


2 


o 

td 



* 


Jt 

QJ 

3fe 





00 

4- 

U- 

CL) 




QJ 





QJ 

o 

to 

-M 

+3 

cr> 

QJ 

CD 

QJ CD 

+3 CD 

+3 CD 

QJ 

CD 

QJ 

CD 

+-> CD 


C 


rd 

cr 

4-> 

C 

4J C 

f o cr 

rd cr 

+3 

cr 

4-> 

sz 

cr 


rd 

!— 

r— 

•r— 

<0 

•r*~ 

fO T- 

i— •*— 

r— •*— 

rd 

•r— 

rd 

•r— 

■r— # r- 

00 

t- 


CU 

>> 

j— 

>> 

«— >> 

cu >> 

CU Js 

r— 

>0 

r— 

>> 

CU 


-M 

+3 


S- 

cu 

C- 

cu l- 

u. 

u. 

Q_ 

u 

CU 

U 

U- 

os 


C 

“O 

(T? 


rd 

rd 

-o rd 

"O QJ 


rd 


rd 

“O ro 

LU 

r— 

QJ 

QJ 

> 

QJ 

> 

QJ > 

QJ > 

QJ > 

QJ 

> 

QJ 

> 

CU > 

I— 

*d 


> 


+3 


-M 

> 

> 

+3 




> 

o 

■o 

QJ 

S- 

r— 

rd 

r— 

rd i— 

V r— 

U. CM 

QJ 

CM 

rd 

CM 

U. CM 

s 

e-r" 

o 

i— 

3 

i— 

r— 

i — 

r- r— 

3 i— 

3 CM 

i— 

CM 

r- - 

CM 

3 CM 


LU 

O 

2T 

LU 

Lu sr 

O ^ 

O '2Z 

LU 


Lu 

2T 

o zz. 

nr 

• 

• 



• 


• 

• 

• 

• 


• 


• 

o 

1 — » 

CM 

i — 


CM 


co 

*3* 

LO 

CO 


r^. 


CO 





QJ 




_£T 




+3 


»— 




rd 

QJ 

4— 


X) 

-C 

O QJ “O 


O 

+3 

xr qj 


C C 

o 

to +3+3 


o 

c 

+3 cu rd 


rd *r- 

rd 

cr l c 


4-> 


QJ +3 O *r- 


QJ rd 

o 

E 4- B 


r«— £ 

-L3 

QJ S- rd 


*r— L. 

C 

* — O to i— 


M- O 

•»— 

<D 4- C 


4— 


o u 


>> L0 

u 

QJ X •«- *r— 


s- c 

o 

r*r+3 a 


rd rd 

+3 

+3 S- rd O 


C S- 

o 

+3 3 U 



QJ 

O) rd CT+3 • 


-Q 

> 

aj E oj o +-> 


r— 


+3 jz cr 


rd rd 

QJ 

rd +3 E +-> CD 


-M 

C 

C 3 U- E 


E c 

o 

<u oj o qj 


O QJ 


C -r- U r- 

LU 

S*- E • 

1 /) • 

QJ U JQ +3 a) 

00 

4- 0J X 

QJ U- 

CD C 

o 

r— ‘ •!— 

U O 

4— p— QJ QJ 

cu 

TMU L 

3 -M 

2T 4~ •»- U +3 

QC 

rd -M 

XJ U 

LU QJ 3 S- rd 

ro 

0 ) S~ rd 

QJ QJ 

CD O cr 3 r— 

Q_ 

QC o E 

QC > 

a: a qj u a 


LU 

sr 

QC 

LU 


M 

h- 

> 


H- 

a 

ro 

zr 

ro 

<c 

o 

LU 

o 

LU 

LU 

CD 

CC 

Of 

QC 

as 


9.20 










£ 


£ 





£ 









O 


o 





o 









*r— 


*r- 





*r— 









P 


4-> 





P 









O 


U 





O 









ai 


CD 





CD 


1— 



oo 


CM 


CL 


CL 



00 

CM 

CL 


00 


r— 

oo 

oo 

CM 

r— 

to 

r— 

to 


r— 

oo 

CM 

to 

CO 

LU 


1 

1 


I 

| 

£ 

1 

£ 


i 

i 

i 

£ 

i 

p 


Lu 

LU 

LL. 

LU 

LU 

►— 4 

Lu 

i — f 


LU 

LU 

LU 

1—4 

Lu 





to 

to 















CD 

CD 









~a 






£ 

£ 








to 

£ 






3 

3 








CD 

O 






P 

P 








£ 

5 





>> 

O 

O 







>> 

3 


CD 





3 

3 








P 

CD 

£ 




c 

£ 

£ 







c 

O 

£ 

O 




o 

P 

-P 







o 

3 

O 






to 

to 








£ 


£ 




to 

-Q 

JO 







to 

P 

£ 

03 




CD 

3 

3 







CD 

to 

03 

-£ 




£ 

to 

to 







£ 

JD 

-£ 

P 




3 









3 

3 

P 




>> 

P 

-Q 

-a 






>> 

P 

to 


CD 




U 

£ 

£ 







O 


CD 

£ 



c 

3 

LU 

LU 






IE 

3 

T3 

£ 

O 



o 

£ 








o 

£ 

£ 

O 

E 




P 

“O 

"O 







P 

LU 

E 



J c 

OJ 

to 

£ 

£ 



£ 



CD 

to 



to 


p 

£ 

-Q 

fd 

<d 



03 



£ 

-Q 

-a 

to 

•i— 


•r— 

3 

3 





r— 



3 

.3 

£ 

«t— 




-P 

to 

P 

P 



3 



P 

to 

05 


>> 



o 


03 

03 

P 


cn 

£ 


U 



>0 

03 


X 

3 

“U 

CD 

CD 

O 


£ 

05 


3 

“O 

P 

03 

£ 

oo 

•r* 

£ 

£ 

a. 

Q. 

> 


•i— 

r— * 


£ 

£ 

cd 

£ 

£ 

o 

£ 

P 

LU 

CD 

CD 

•r~ 


CO 

3 


P 

LU 

CD 

£ 

03 

»— 1 

P 

</> 


OH 

CC 

CU 

p 

1 

cn 


CO 


CL 

03 


1 — 

03 

_Q 

*o 




o 

£ 

£ 


JD 

-o 

CD 


cn 

00 

£ 

3 

£ 

r— 

CM 

o 

> 

o 

•r- 


3 

£ 

C£ 

CO 

£ 

l-H 


to 

05 



£ 

*r- 

£ 

to 

JC 

to 

03 


£ 

•i — 

cn 

U 



•» 


CD 

CL 



P 




•r— 

> 

LU 


P 

p 

-P 

-P 

N 


to 

to 

•f— 

p 

P 

p 

P 

•r— 

H 

£ 

£ 

£ 

£ 

£ 

1 

o 

•f— 

•r— 

2 

£ 

£ 

£ 

£ 

CD “O 

CJ> 

p 

o3 

1 0 

03 

fd 

£ 

£ 




03 

rd 

03 

03 

U £ 

<c 

CD 

-P 

p 

-P 

*P 

o 

CD 

p 

P 

r— 

P 

p 

P 

p 

CD O 

c2 

E 

00 

00 

OO 

00 

z 

CM 

X 

X 

CD 

OO 

00 

OO 

oo 

OC 2 

<c 

E 









£ 






nz 

>> 

• 

• 

• 

• 

• 

• 

• 

• 

rd 

• 

• 

.• 

• 

• 

o 

OO 

i — 

00 

CO 


r— 

CM 

p— 

CM 

Q_ 

i — 

CM 

00 

i — 

CM 






X 




£ 


•r— 




o 


£ 




•r- 


P 




P 


03 O 



P (O 

CJ 

CO 

E p p 



£ 03 

3 


-Q £ O 

03 


03 _£ 

"O 

f( 

3 *f- 

£ 


£ 

CD 


to CD 

P CD 


•i- _£ T3 

£ 

<c 

£ £ • 

03 j£ 


E o CD 



P 3 X 

03 P 


£ *r~ P 

£ 

• 

£ 03 P *r- 

cn o 


CD JT 03 

03 


<D a O £ 

£ O £ 


P S CD 

•r* 

x 

E 3 P 

•«— P fd 


CD CL 

to • 


CD CD £ 03 

p 


■o X <D 

to X 

E 

*— X? P E 

£ X £ 


•i- £ 

3 •«- 

CD 

CD O CO 

fd fd *p*“ 


CD £ 

03 £ 

r“— 

£ -Q to 

P £ 


-£ P “O 

O P 

-Q 

'd" 3 to 

to £ £ 


P 03 £ 

03 

O 

X 03 CO CD 

03 O 


E 03 

03 E 

£ 

^ £ 

tO T— 


to 


Cl 

-£ CD P 

£ CD P 


cd u cn 

to U 


03 P P P 

CD £ *.— 


P -1“ £ • 

E **“* 

X 

•r— fd *r— 

P o to 


3 £ •!— CD 

£ £ 

•r* 

tO ^ *r- P 

U O 


CLP CL £ 

O P 

£ 

CD £ tO 

fd £ CL 


E CD CL 3 

P CD 

P 

cn-o cl 

£ -r- 


O E 03 P 


03 

£ CD O r— 

03 P 


U E »— U 

CD E 

E 

CD P £ CD 

-£ £ £ 

LU 

>> £ 3 

CL X 


E 03 CL C 

U O CD 

OO 

h— to CD £ 

to 

to • 

•r— CL 03 

•«- £ • 

o 

a > p 

r 

CD X 

LU U 03 Q- 

CO P CD X 

CL 

X 03 O to 

»-h 03 

> 

OC O 

CD -t- P 03 

CC 

_U -Q 

P 

r— £ 

O to CD CD 

> (O P £ 

ZD 

CQ P £ 3 

LU £ 

o o 

h wrr 

O O * £ 

Q_ 

OO O 03 tO 

CO O 

co P 

CO 03 P P 

s: CLT3 03 


LU 

►—4 

P 

Q 

3C 

LU 

LU 

> 

i 

P 

X 

t— i 

>• 


ID 

_J 


-J 

o 

QC 

O 

CQ 

LU 

o 

1 — 

! — 

C£ 

CO 

CO 

oo 

oo 

OO 


9.21 



ROUTINE PURPOSE CHARACTERISTICS TEST 


o 


-M 

0 

01 

CL r— VO LO 

tO i— ^ n r» n CO 

E 8 I I S 8 i 

9— t LL. U„ II U- Ll Ll 


CNJ CO LO VO CO CPV t— 

I i i i 1 I i i 

U_ U_ U-U.Ll.U- Ll> Uu 







e 









o 

e 








•r* 

o 








+-> 

•r— 


sz 






CL 

-M 

-E 

to 






O 

CL 

00 

to 







O 

fd 

I— 






s- 


r— 

to 






o 

cn 

■CO 


xj 



>> 


•M 

e 



03 





O 



o 

L- 



E 

>> 

03 

JC 

o 


03 

to 


o 


> 

u 


-M 

•M 

-M 

to 


e 

E 


4-> 

-X 

e 

O 

«M 

e 

o 

01 

<d 

4— 

cn 

3 

> 

a; 

o 


o> 

a> 

CU 

•r— 

o 

•r— 

CO 

•r» 

e 

•r— 

to 

r— 

L 

u 

CL 


-M 

o 

Ui 




e 


to 

CL 

•r— 


•M 

E 

C 

03 

O 

4-> 

O 

-M 

-a 

O 

o 

o 


JL 

H3 


CL 

e 

o 

•r— 

•I— 

«M 

cu 

-a 

-x 

a 

rd 


*M 

*M 

O 

N 


u 




CL 

CL 

> 


01 

a» 

03 

o> 

CU 

o 

o 

•r— 

3: 

r*-» 

-E 

e 

e 

r— 



CL 

03 

a. 

u 

*r- 

•r* 

CL 

-M 

-M 


E 

•*— 


r_ 

r— 

•r* 

to 

03 

O 


-M 

fd 

-X 

-X 

-M 

CU 

CU 

L- 

r— 

r— 

4-> 

U 

a 

i— 

CL 

CL 

0) 

r— 

3 

fd 

3 

3 

3 

CU 

CU 

rM 

<c 


o 

CO 

CO 


cu 

CU 

• 

* 

CM 

.# 

• 

Osl 

• 

oo 

• 

• 

u> 

• 

• 

CM 






to 




w* 



-M 




to 



-a o 




>0 



U 3 




fd 



rd S- 







O *M 


03 


*•- H 



to 


to -E 


O 03 0) 



E 


to -M 


> E to 



<u o 


03 tO 



o> 


+J CJ 


E O 4-> 

to 

O -MO 

E 


E 


4— «M C 



•i— 


•1- XJ 


4- 03 

u 

-M £? s- -M 

to 


O E 


•r W E 

o 

3 rd O C 

to 


CL rd 


+J ■r— 03 

4- 

U E w O 

o 




lO X T— 


CJ .r w U 

s- 


03 to 


rd 03 « 

V- • 

X S- 0) 

o 


j— T3 


r- 03 

CU X 

cu cl u -a 



X) L- 

• 

fd r— 03 “O 

> -r- 

O E -M 


fd rd 

to 

-a fd -m o 

r- L. 

a; to s- fd 

o 


-M O 

L 

O U fd E 

O -M 

r l> a 

o 



03 

E O f— 

to <d 

-M »r- E 

i- 


CU CQ -M 

r- ac 

E 

r— O 



-E 

E 

to O 

E 

CO XI fd •»— 



+-> X3 

•*— 

E ai 

O X) 

r- O) 5- -M 

«— i 


E 

O 

C £= o i- 

•r- E 

O 3 -M rd 

+ 


to rd 

CL 

O *M 4— 03 

-M fd 

S- O C E 

c 


-a 


4- «M 

fd X) 

-M £ 03 E 

' — 


fd Cu 4- 

to E to C 

3 

C£ UO 

tr 


03 

O 

E O *r- •«— 

O' CJ 

O 4-> 4— 


• 

L. a> 


rd S- X 1 

cu T- 

o CO E 


to 

E 

>> 

i- 4- rd E 

l. 

a. 

to 

a> 

h- *r- 

rd 

-M O 

L. 4J 

<(0*0 

3-C 

Z 3 

s- 

X r— E 

rd 03 

^ «C •*- cu 

O 

o 

1— 1 o 

s- 

Ll. •«— rd 

CU E 

<n _j > o> 

•f— 

s- 

o »— 

rd 

2: S- X3 -E 

E E 

COO O 

s- 

<d 

CL 


<4J O+J 

•r- >» 

t. to 

rd 

cu 

CQ O 

E 

tt d i-*r- 

— J to 

CO s CL 3 

> 

to 

h- 4- 

rd 

h E cnj 



9.22 


TRANSF TRANSF transfers specified blocks 
of rows and columns from one array 
to another array. 



E 

O 


1 — 


O 

O 

O 

CO 


r— 

r— 


LU 


1 

% 

I 

h- 


U. 

Ll_ 

U- 








CD 

P 




E 

O 




CO 





CL 

P 





<0 




P 

CD 




O 

CL 



<u 


<U 



x: 

CD 

E 



p 

L 





3 

P 



E 

P 

CO 

i— - 


•r— 

U 

CO 

CD 



3 

r- 

E 


f— 

E 


CO 


01 

P 

o 

CL 


e 

CO 

p 



co 

-Q 


P 


CL 

3 

CD 

O 



CO 

e 



-o 


3 

CD 


CP 

-a 

P 

E 


> 

e 

U 

3 


e 

o 

3 

4-> 

LD 

3 

u 

E 

U 

O 

U 

<p 

P 

3 

i— i 


CO 

CO 

E 

1 — 

CO 


-Q 

P 

oo 


e 

3 

CO 

1— 1 

e 

o 

V) 

-Q 

o: 

•p— 



3 

LU 


P 

-O i— 

CO 

I— 

P 

CO 

E CD 


o 

E 

e 

•i— E 

-o 

c 

CP 

•r— 

-E <0 

E 

cx 

E 

U_ 

H- CL 

LU 

<c 

CD 




□c 

r— 



* 

E> 

LU 

CO 

JD 

o 


n-Cvjr0^tLO^D00O> 
118 18 118 
(jl.Ll.U-U-Ll.Li_U.Li- 


o 

CP 

D- 

CO 


i/> 


O 



•r— 


<D 


P 


-E 


U E E 

in 

P 


CD *i- *r- 

•r* 

CD 


f— CD CD r— 

X 

E 1 -E 


■4- XT XT CO 

E CO 

O P H- 4-> 


CD P P X> 

O 

EX: CD 


-o o 

•r- r— 

P CD • CO 


E 0«— 

4-> CO 

•p— CD 


P O P 0> 

O -Q 

*3 E CD 


E E 

CD O 

8 T3 E 


0) E "O 

r— r— 

h di E CO 


E O CD 

4- O) 

x: o 


a; i— E 

CD 

CO P i CD 


r- ©<P E 

3 0) • 

p 


CD E •«— O 

jz co 

•r- CO CD P 


fO CO P 

r— P -r- 

x ax o 


CD Q- O CO 

C0 X 

U P 


P CL E 

-a E <o 

CO co 00 


CO -O CO 

o o 

P r- O P 


i— CD r— E 

E E »— 

U CLP -r- 


a> o-p 

4- CO 

CO EX) 


E X) E 

co a 

E*0*r 


E 3 O 3 • 

EP O 

pe a> 


E O r— E 

E E »— 

X flO T3 E • 

z: 

O OOr- O 

O CD 

CD CD -r- O 

<c 

4- CD CO-»- 

4- E CD 

C P E E 

QC 

W r r- P 

w or 

O </>•»— CD 

CD 

E P fO P*r- 

E CLP 

C -o 3 CO N 

O 

CO E •?- CO 

<0 E 

Q_ E £ 

CX 

E E -r- E O 

E 5 o 

ZOO 0)0 

Q- 

J— t- 4- -r- CL 

J— O P 

Z3 % CO E P 


LU 


s 


CO 

«k 

o 

o *» 

c 1—5 
CL. •* 

3E 


9.23 



9.24 


fv* 

CM 

I 


CO 

oo 

1 

Ll 


oo 


oo 


oo 


co 

o 

Cl. 

CSC 


3 

<d 

o 


<0 

o 


i- 

u 


a» 

JZ 


i- 

o 

4- 


3 

0 

_a 

s- 

0 > 

CL 

CL 

3 

CD 

oG 

+-> 

to 

01 
+3 
id 


id 

a 


2 

Q- 


3 ) 

§ 


QC 

a. 

CL 


co 

z 

o 


Z Q CD 

M CQ •— >4 


CD 

•— i 

Ll 

z 

o 

o 


CO 

4 -> 


^ g 

•p" (D 

^ OJ 
CO 


4 -> 

* 

CD 

4 J 

ro 


S- 

td 


ill 


*— a) 
-Q -Q 

2 ~ 
a. <a 


CO 

xT 


00 


OO 

00 


CO 


p— 

CM 

OO 


un 

to 


OQ 

oo 

CM 

OO 

oo 


CM 


f— 

p— ^ 

r— 

f— 

p— 

f— 

r— 

p— - 

s 

g 

8 

1 


r 


8 

1 

i 

0 

i 

i 

1 

0 

LL 

lL 

Ll- 

Ll 


Ll 


Li- 

U- 

Ll 

Ll 

ll 

Ll 

Ll 

Ll 













X? 

-o 














0 ) 

CD 














-L> 

4-> 














O 

U 













XJ 


01 













CD 

4- 

4- 

44 












44 

4- 

4- 

O 












O 

«£ 

<d 

id 

G 












CO 

-L> 

to 












ci- 

•p— 

O 

•*— 












io 

4-> 

C 












"O 

to 

C 

to 

c 











CD 

•r- 

CD 

°r— 

aj 













E 


E 











u 

*L> 

CD 

•M 

CD 









-o 

~a 

CD 

3 

r— 

C 

p— 









<v 

CD 

4- 

CD 

CD 

Q> 

CD 









+J 

-L> 

4~ 

E 


E 










u 

U 

<0 

CD 

CD 

<D 

CD 








3 

a> 

CD 


i— 

-M 


4-> 








<D 

4- 

4— 

4-> 

0 ) 

id 


id 








-P 

4- 

4- 

O 


r— 


r— 








U 

id 

<d 

C 

0 ) 

CL 

0 ) 

CL 








CD 




4~> 











4- 

co 

-M 

to 

id 

+-> 

id 









4— 

*f— 

o 


r— 

to 

r— 

to 








id 


c 


CL 

S- 

a. 

■G 

CO 








oo 


oo 


•r- 


•p— 

JC 







to 


to 



Ll 

4-> 

Ll 

-P 


co 





*r— 

4-> 

•r- 

-L> 

to 


to 


3 


x: 






c 


c 

S- 

«k 

L. 


•r— 







OO 

<d 

oo 

fl} 

•r- 

XJ 

•f— 

-o 

* 


-a 






E 


E 

Ll 

CD 

Ll 

CD 


•r- 





-P 

CD 

*L> 

a> 


4-> 


-L> 

-P 


3L 





C 

r— 

C 

r— 


U 

•k 

U 

C 







CD 

Si 

CD 

X5 

a> 

3 

CD 

0 ) 


-P 





E 


E 


CD 

4- 

0 ) 

4— 

c. 


C 





CD 

CD 

(D 

CD 

4-> 

4- 

-L> 

4- 

CD 


cd 





r— 

+-> 

i— — 

4-> 

U 

id 

a 

10 

4- 


c. 





CD 

Id 

CD 

id 

CD 


3 


4- 


a> 






r— 

r— • 

4- 

+-> 

4- 

-L> 

•r~ 


4- 





CD 

CL 

CD 

CL 

4- 

o 

4- 

O 

3 


4— 





4-J 


-M 

io 

c 

10 

G 



«r— 





id 

4-> 

id 

-M 




©3 


3 





i— 

CO 

p— 

to 

to 

to 

to 

to 








OL 

S- 

CL 

S- 

•r— 

«r— 

•f— 

•r- 

CD 


od 






•P-. 


•r- 





E 






-L> 

ll 

+-> 

Ll 

3 

3 

*3 

3 

id 


CD 





to 


to 






10 


E 





s- 

•> 

S- 


+-> 

-D 

4-> 

44 


*d 





•r- 

-o 

•r— 

-a 

C 

C 

c 

C 

JC 


CO 





U- 

CD 

Ll 

CD 

fl) 

CD 

m 

CD 

•p 








+J 


4J 

F 

E 

p 

E 

•1— 

g 

JC 






U 

ct 

u 

CD 

CD 

0) 

CD 


o 

•P 






<u 

-o 


r— 

r— 

r— 

r— 

•*— 

*f— 

CO 




CD 

4- 

<D 

4- 

CD 

CD 

CD 

CD 


•p 

3: 

to 





4- 

4-> 

4- 






•r- 


CD 




U 

id 

U 

id 

01 

CD 

3 

CD 

c 

3 


c 




CD 


3 


4-> 

4-> 

-M 

-M 

o 

c 


4- 




4— 

-D 

4- 


id 

Id 

id 

id 

*r— 

o 

CO 

4— 




ci- 

O 

4- 

o 

r— 

i— » 

r— 

i— 

-P 

u 

CO 





id 

C 

id 

c 

CL 

CL 

Q- 

CL 

•r* 


tt) 

*P 












3 

>> 

c 

to 




to 

to 

to 

to 

m 


<k 


c 

u 

4- 





°f— 

•r— 

*r- 

•I— 

o 

o 

o 

o 

o 

Id 

4- 

+-> 












u 

3 

•r* 

C 




•3 

X> 

*3 

•3 

0 L 


TL. 

-1L. 


c 

4-> 

•r— 



i — 

i 








>> 

3 

CO 

G 

-C 


CM 

+-> 

+-> 

+J 

4-> 

O 

O 

O 

O 

L 

o 


CL 

-P 


CM 

C 

c 

C 

c 

•*— 


•r— 

•r- 

<o 

X> 

VJ> 

to 

•I— 


'SC 

ft 

Q) 

CD 

CD 

4-> 

X> 

+-> 

-L> 

3 


c 


3c 



E 

E 

E 

E 

id 

id 

id 

id 

c 

-P 

•l— 

-P 




CD 

CD 

CD 

CD 

s- 

s- 

S-. 

u. 

3 

c 

L 

C 

to 

o> 

CO 

r* 

r- 

p- 

r— 




3 

o 

CD 

a. 

CD 

E 

c 

c 

CD 

Oi 

CD 

a) 

p— 

1 — 

r— 

cu 

J3 

SL 

CO 

S- 

<D 

• 1 — 

•t— 





id 

id 

id 

Id -L> 


<u 


CD 

r— 

>> 

>> 

<D 

CD 

CD 

CD 

•r— 

■r* 

•«— 

•r- u 

V 

4- 

o> 

4 - 

XI 

s- 

sJ 

-M 

4-> 

4 > 

4-> 

X 

X 

X 

X CD 

E 

4- 

E 

4- 

o 

<d 

id 

id 

id 

id 

id 

id 

id 

id 

Id 4- 

m 

•i— 

id 

•r* 

s- 

> 

> 

r— 

r— 

r— 

r— 

•r— 

•r— 

•r 

•p— ^ 

tO 

3 

-40 

3 

CL 

CM 


CL 

CL 

CL 

CL 

CQ 

CO 

CQ 

CQ id 



<*«— S 


CD 

CM 

r— 

s 





#»— s 



r— 

CM 

oo 

«Sf 

•P 

Z 

z 

r— 

CM 

OO 


in 

to 


oo 


<o 

CL 

« 

CM 


CD -Q 


9.24 



ROUTINE PURPOSE CHARACTERISTICS TEST 
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10.0 


PROGRAM VALIDATION 


The program s BUCLASP2, was validated by running various data sets such 
that all major logic paths were tested. The validation is divided 

into two categories: 

a. Functional Tests 

b. Inspection 

Listed below Is the functional tests and their characteristics. 

Section 9.3 states for each program/subroutl ne the test(s) used. 


10.1 



10.1 


Input Option Functions] Test Correlation 


Field 

Option 

Test 

1 

Data Set Title 

F-l 

1 

IPC(l) 

3 0 

F-l 


IPC(l) 

* 1 

F-81 

2 

IPC(2) 

* 0 

F-l 


IPC(2) 

- 1 

F-82 

3 

IPC(3) 

3 0 

F-l 


IPC(3) 

3 1 

F-83 

4 

IPC(4) 

* 0 

F-l 


I PC (4) 

3 1 

F-83 

5 

IPC(5) 

3 0 

F-T 


IPC(5) 

3 1 

F-84 

6 

IPC(6) 

3 0 

F-l 


IPC(6) 

3 1 

F-85 

7 

I PC (7.) 

3 0 

F-l 


IPC(7) 

3 1 

F-86 

8 

IPC{8) 

3 0 

F-l 


I PC (8) 

3 1 

F-87 

1 

JPC(l) 

3 0 

F-l 


JPC(l) 

3 1 

F-41 

2 

JPC(2) 

* 0 

F-36 


JPC(2) 

3 1 

F-l 

3 

JPC(3) 

3 0 

F-l 


JPC(3) 

3 1 

F-41 

4 

JPC(4) 

3 n 1 <_ n < 100 

F-89 

5 

JPC{5) 

3 k 4 < k < 10 

F-37 

1 

Number of Nodes 

F-l 

2 

Number of Beams 

F-41 

3 

Number of Flat Plates 

F-41 

4 

Number of Curved Plates 

F-41 


10.2 
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Continued 


• • 9 


Card 

Field 

Option 

Test 

4 Cont 

« • 4 




5 

Curved Panel Option 




« 0 

F-1 



« 1 

F-10 


6 

Curved Panel Angle 

F-10 


7 

Section Length 

F-1 


8 

Load or Strain Value 

F-1 9 


9 

Biaxial Ratio 

F-1 5 

5 

1 

Load Option 




« 1 

F-ll 



« 2 

F-28 



- 3 

F-19 


2 

Wave Number Option 




- 1 

F-37 



- 2 

F-1 



* 3 

F-36 



- 4 

F-39 


3 

Loop Start Value 

F-1 


4 

Loop End Value 

F-1 


5 

Lower Limit for Root Search 

F-35 


6 

Upper Limit for Root Search 

F-35 


7 

Default Number of Subdivision 




■ blank 

F-1 



* n 

F-37 


8 

Lower Bound Load 

F-37 


9 

Upper Bound Load 

F-37 

6 

1-16 

H-list Values 

F-36 

7 

1-14 

H-list Values 

F-36 

T-l 

variable 

Thickness Table 

F-37 
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10.1 


Continued 


.* ® ® 


Card 

Field 

Option 

Test 

T-2 

variable 

Material Table 

F-37 

C 

1-9 


F-27 

E 

1-9 


F-27 

N 

1 

Letter N 

F-2 


2 

User Node Number 

F-2 


3 

Nodal Boundary Condition 




* blank 

F-2 



« 1 

F-2 



- 2 

F-l 



* 3 

F-l 



* 4 

F-l 1 


4 

Y-coordinate 

F-2 


5 

Z-coordinate 

F-2 

D 

1 

Letter "D" 

F-l 


2 

Number of Substructures 

F-1 


3-8 

Substructure Node Pairs 




Start Only 

F-l 



Start and End Only 

F-38 



Start, Repeat and End 

F-22 

S 

1 

Letter "S" 

F-ll 


2 

User Node ID 

F-l 1 


3 

W Component 

F-ll 


4 

0 Component 

F-ll 


5 

v Component 

F-ll 


6 

u Component 

F-ll 

F 

1 

Letter "F" 

F-22 


2 

User Plate ID 

F-22 
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Conti nued 


® » « 


Card Field 

Option 

Test 

F (Cont .. .) 



3 

0FF1 

F-22 

4 

0FF2 

F-22 

5 

0FF3 

F-22 

6 

0FF4 

F-22 

P 1 

Letter P 

F-l 

2 

User Plate ID 

F-l 

3 

Node I 

F-l 

4 

Node J 

F-l 

5 

Plate Type 



- 1 Flat 

F-l 


= 2 Curved 

F-30 

6 

Element Offset 



* 0 No 

F-l 


* 1 Yes 

F-22 

7 

Number of Laminas 

F-l 

8 

Input Mode 



* 0 P-1, or P-2 card 

F-l 


* 1 Table 

F-37 

9 

N ^2 Load Switch 

F-38 


. « 0 affected 

F-38 


« 1 not affected 

F-36 

10 

Number of Subdivisions 

F-37 

11 

Type of Material Properties 



* 0 Engineering Constants 

F-37 


* 1 Q matrix 

F-37 

P-1 1 

T - 

F-l 

2 

E 11 

F-l 

3 

E 22 

F-l 

4 

RNUA 

F-l 

5 

G 12 

F-l 
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Card 

Field 

Option 

Test 

P-2 

1 

T 

F-37 


2 

Oil 

F-37 


3 


F-37 


4 

q 22 

F-37 


5 

«66 

F-37 

B 

1 

Letter "B" 

F-23 


2 

User Beam ID 

F-23 


3 

Beam Node 

F-23 


4 

Type of Beam 




= 1 

F-25 



* 2 

F-24 



* 3 

F-26 


5 

Number of Layers 

F-23 


6 

Input Mode 




= 0 B-l » B-2 or B-3 card 

F-23 



= 1 By tables 

F-39 


7 

Beam Angle 

F-23 


8 

Beam Area 

F-23 

B-l 

1 

E-modulus 

F-25 


2 

G-modulus 

F-25 


3 

Moment About Local y-axis 

F-25 


4 

Moment About Local z-axis 

F-25 


5 

Warping Constant 

F-25 


6 

Torsion Constant 

F-25 


7 

Y offset 

F-25 


8 

Z offset 

F-25 

B-2 

1 

Thickness 

F-23 


2 

E-modulus 

F-23 


3 

G-modulus 

F-23 


4 

Beam Width 

F-23 
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Continued 


• ® 9 


Card 

Field 

Option 

Test 

B-3 

1 

Radius 

F-26 


2 

E -modulus 

F-26 


3 

G-modulus 

F-26 

L 

1 

Letter "L" 

F-l 
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Functional Test Descriptions 


Test 


Characteristics 

F-l 

1. 

Flat Panel at 90° SS - Free 


2. 

Buckling and eigenvector solution 


3. 

M0PT = 2 


4. 

Flat Plates with 



5. 

a. left node reduced out 

b. both nodes interior 

c. right node reduced out 
Input by P-1 cards 



6. 

Simply-supported node 


F-2 

1 . 

Flat Panel at 90° C - Free 



2. 

Clamped node 



3. 

Free node 


F-3 

1 . 

Flat Panel at 0° angle SS - Free 


F-4 

1 . 

Flat Panel at 0° angle C - Free 


F-b 

1 . 

Panel at angle (non-orthogonal) SS - 

Free B.C. 

F-6 

1 . 

Panel at angle (non-orthogonal ) C - 

Free B.C. 

F-7 

1 . 

Flat plate with non-zero B matrix 


F-8 

1 . 

Curved panel at angle (1) 


F-9 

1 . 

Curved panel at angle (2) 


F.-10 

1 . 

Repeat Block Curved Panel 


F-l 1 

1 . 

Nj 1 varying 



2. 

N 22 i Input 
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F-ll 


F-12 


F-13 


F-14 


F-15 


F-16 


F-17 


Conti nued 


3. First Plate Affected 

4. J th Plate Affected 

5. Sprung 6.C. 

(J is the critical element for panel upper bound.) 


1. Njj Varying 

2. N 22i Input 

3. First Plate Affected 

4. J th Plate Not Affected 

1. N-j-| Varying 

2. N 22i Input 

3. First Plate Not AFfected 

4. J th Plate Affected 

1. N-ji Varying 

2. N 22i Input 

3. First Plate Not Affected 

4. J th Plate Not Affected 


1. Biaxial Ratio 

2. First Plate AFfected 

3. J th Plate Affected 


1. Biaxial Ratio 

2. First Plate Affected 

3. J th Plate Not Affected 

1. Biaxial Ratio 

2. First Plate Not Affected 

3. J th Plate Affected 

10.9 
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Continued 




Test 

F-18 


F-19 


F-20 


F-21 

F-22 

F-23 

F-24 

F-25 

F-26 


Characteristics 

1. Biaxial Ratio 

2. First Plate Not Affected 

3. 0 th Plate Not Affected 

1. N 22 Varying 

2. EPS 11 Input 

3. First Plate Affected 

4. J th Plate Affected 

1. N 22 Varying 

2. EPS 11 Input 

3. First Plate Not Affected 

4. J th Plate Affected 

1. Mix of Curved and Flat Plates 

1. Offsets on Plates 

2. Block Interrelationship Node Pairs 

3. Plates with Different A, B & D Matrices 

1 . Rectangular Beam with 2 Layers 

2. Angle of Beam is Non-Trivlal 

1 . Rectangular Beam with 1 Lamina 
1 . General Beam Element 
1. Circular Beam with 3 Lamlnas 
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Continued 
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Test 
F -27 


F-28 


F-29 


F-30 


F-31 


F-32 


F-33 


Characteristics 

1 . Plate with 3 Layers 

2. Rectangular Beam with 15 Layers 

3. Similar Plates with Different Boundary Conditions 

4. Similar Plates with Same Boundary Conditions 
with Different Widths 

5. Non-Neglected Plate Element with Non-Interior 
Node 

1. Flat Plate with Zero B-Matrix 

2. N 22 Varying 

3. Non-Zero N^ Side Load Input 

1. Flat Plate with Zero B-Matrix 

2. N^ Varying 

1. Curved Plate with Non-Zero B-Matrix 

2. Non-Zero N^ Side Load Input 

3. N 22 Varying 

1. Curved Plate with Non-Zero B-Matrix 

2. Non-Zero N 22 Side Load Input 

3. Nii Varying 

1 . Determinant Always Positive 

2. 0/2 Buckling Condition 

1. Curved Elements with Zero B-Matrix and Different 
Radii 

2. Similar Plates with Same Boundary Conditions 
with Same Widths 

3. Similar Plates with Same Boundary Conditions with 
Same Widths and Same Spring Sti ffness 
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Continued 
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Test 

F-34 

F-35 

F-36 


F-37 


F-38 


F-39 


F-40 


F-41 


Characteristics 

1 . Curved Plate with Zero B-Matrlx, N 22 

Varying 

1. Multiple Root Searching 

1. H 22 Varying, Input Case with Some 

Unaffected Elements. M0PT * 3 Used and 
Buckling Only Option Used 

1. Flat Plate with Blj f 0 with Input Upper Bound 
and M0PT * 1. Qij of One Plate Input with P-2 
Card, Qij of Another Plate Input with Table and 
Card. E, v Input by Table, JPC(5) « 6; 

Default Subdivisions * 10; One Element Set * 5 

1. Two Identical Plates with Different Spring 
Stiffnesses at the Same Sprung Sides Going 
Through the Upper Bound Calculations. Only 
Start and End Substructures are Used for This 
Problem. 

1. One Beam is Circular with One Layer, the Other 
is a General Beam with Offsets. The Circular 
Beam is Input with Table. M0PT * 4 is Used. 

1. Two Laminated Plates with Same Loads and Radii, 
But with Different Widths. These Plates are 
Not Considered Identical From the Stiffness 
Standpoint. 

1 . Too Many Nodes 
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Continued .. . 


Test 

Characteristics 

F-42 

1. 

Too Many Elements 

F-43 

1 . 

Too Many Beams 

F-44 

1 . 

Too Many Plates 

F-4b 

1 . 

2. 

Bad Equivalent Node Data 

Equivalent Nodes with Different Boundary 

F-46 

1 . 

Bad Interrelation Node Data 

F-47 

1 . 

Too Many Values in M-List 

F-48 

1 . 

No L Card 

F-49 

1. 

Non-Element Data Out of Order 

F-bO 

1 . 

Too Many Nodal Spring Stiffnesses 

F-bl 

1 . 

Too Many Plate Offsets 

F-b2 

1 . 

Too Many Interrelationship Node Pairs 

F-b3 

1 . 

Too Many Equivalent Node Pairs 

F-b4 

1 , 

Bad Table Input 

F-b5 

1 . 

Incorrect Number of Nodes Input 

F-b6 

1 . 

Bad Thickness Table 

F-S7 

1 . 

Too Many P Cards 

F-b8 

1 . 

Too Many B Cards Input 

F-59 

1 . 

Invalid Load Option 


10.13 
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Continued 


Test 

F-60 

F-bl 

F-62 

F-63 

F-64 

F-bb 

F-6b 

F-b7 

F-68 

F-b9 
F-70 
F-71 
F-72 
F- 73 
F- 74 
F-7S 


Characteristics 

] . Invalid Wave Number Search Option 

1 . Too Many N Cards 

1. Illegal Boundary Condition 

1. Wrong Number of Spring Stiffnesses Input 

1. Bad User ID on Nodal Stiffness Set 

1. Bad User ID in Interrelationship Data 

1. Bad User ID in Equivalent Node Data 

1. Bad User ID in Substructure Definition 

1. Incorrect Number of Substructure Definition 
Node Pairs Input 

1. Substructure Definition Node Pair Out of Order 

1, Interrelationship Between Blocks Invalid 

1. Invalid Node Used on Plate 

1. Invalid Node on Beam Definition 

1. Incorrect Number of Plates Input 

1 . Incorrect Number of Beams Input 

1 . Node in Start Substructure with Attachment to 
Element to Repeat, But Not Part of a Connection 
Node Pair. 

10.14 


10.14 
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Continued 


Test 

F-76 

F-77 

F-78 
F- 79 
F-80 

F-81 
F-82 
F-83 
F-84 
F-85 
F-86 
F -87 
F-88 
F-89 
F-90 


Characteristics 

1 . Block with No Interior Node 

1. Panel with 2 Substructure, But No End Sub 
structure 

1, Boundary Condition on Beam Node 

1. Beam Node Which is Not Equivalenced 

1 . Element with Nodes in Start and End Sub- 
structure 

1. IPC(l) = 1 

1. IPC(2) * 1 

1. IPC(3) * 1 

1. IPC(4) - 1 

1. IPC(5) * 1 

1. IPC(6) » 1 

1. IPC{7) = 1 

1. IPC(8) = 1 

1. JPC(4) - n < 100 and n t 0 

1. Duplicate User ID on Nodes 

2. Duplicate User ID on Beams 

3. Duplicate User ID on Plates 


10.15 
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